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Abstract. Consider a TV X n random matrix Y n = 0"y) where the entries are given by 

ij ^ H > 

the X™j being centered, independent and identically distributed random variables with 
unit variance and (cr^n); 1 < i < N, 1 < j < n) being an array of numbers we shall refer 
to as a variance profile. We study in this article the fluctuations of the random variable 

logdet (Y n Y* + pi N ) 

where Y* is the Hermitian adjoint of Y and p > is an additional parameter. We prove 
that when centered and properly rcscaled, this random variable satisfies a Central Limit 
Theorem (CLT) and has a Gaussian limit whose parameters are identified. A complete 
description of the scaling parameter is given; in particular it is shown that an additional 
term appears in this parameter in the case where the 4 th moment of the JTjj's differs 
from the 4 th moment of a Gaussian random variable. Such a CLT is of interest in the 
field of wireless communications. 
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1. Introduction 

The model and the statistics. Consider a N x n random matrix Y n = (YQ) whose entries 
are given by 

Y- = ^lx- , (1.1) 

where (<Tij(n), 1 < i < N, 1 < j < n) is a uniformly bounded sequence of real numbers, 
and the random variables are complex, centered, independent and identically distributed 
(i.i.d.) with unit variance and finite 8 th moment. Consider the following linear statistics of 
the eigenvalues: 

1 1 N 

l n (p) = - logdet (Y n Y: +pIn) = j^Y. + P) 

i=i 

where In is the N x N identity matrix, p > is a given parameter and the A^'s are the 
eigenvalues of matrix Y n Y* . This functional known as the mutual information for multiple 
antenna radio channels is very popular in wireless communication. Understanding its fluctu- 
ations and in particular being able to approximate its standard deviation is of major interest 
for various purposes such as for instance the computation of the so-called outage probability. 
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Presentation of the results. The purpose of this article is to establish a Central Limit 
Theorem (CLT) for T n {p) whenever n — > oo and ^->c£ (0,oo). 

The centering procedure. It has been proved in Hachem et al. [17j that there exists a sequence 
of deterministic probability measures (7r n ) such that the mathematical expectation ET n (p) 
satisfies: 

EJ n (p) - / log(A + p)n n {dX) ► . 

Moreover, Jlog(A + p)ir n (dX) has a closed form formula (see Section [2T3]) and is easier to 
comput^l than EI„ (whose evaluation would rely on massive Monte-Carlo simulations). For 
these reasons, we study in this article the fluctuations of 

i logdet(F„y„* + pI N ) - J log(p + t)n n (dt) , 

and prove that this quantity properly rescaled converges in distribution toward a Gaussian 
random variable. Although phrased differently, such a centering procedure relying on a 
deterministic equivalent is used in [Tj and [3J- 

In order to prove the CLT, we study separately the quantity N(X n (p) — EX„(p)) from 
which the fluctuations arise and the quantity N(WE n (p) — J log(A + p)tt 71 ( d\)) which yields 
a bias. 



The fluctuations. We shall prove in this paper that the variance 0^ of N(T n (p) — EZ n (p)) 
takes a remarkably simple closed-form expression. In fact, there exists a n x n deterministic 
matrix A n (described in Theorem 13. ip whose entries depend on the variance profile (cr^ ) 
such that the variance takes the form: 

Q 2 n = logdet(I„ -An) + nTrA n , 

where k = E\Xu\ 4 — 2 in the fourth cumulant of the complex variable X\\ and the CLT 
expresses as: 

N c 
— (I„-EX„) -±_+^(0,l). 

H n n—too 

In the case where k — (which happens if Xij is a complex gaussian random variable for 
instance), the variance has the log-form 8^ = logdet(/„ — A n ). This has already been 
noticed for different models in the engineering literature by Moustakas et al. [22j . Taricco 
[29] , See also Hachem et al. in [14]. 

The bias. It is proved in this paper that there exists a deterministic quantity B n (p) (described 
in Theorem 13.31) such that: 



N EI n (p) - / log(A + p)n n ( dX) - B n {p) ► 

V J J n ^°° 

If k — 0, then B n (p) — and there is no bias in the CLT. 



^especially in the important case where the variance profile is separable, i.e., where crfj{n) is written as 
<Tij(n) = di(n)dj(n). 
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About the literature. Central limit theorems have been widely studied for various models 
of random matrices and for various classes of linear statistics of the eigenvalues in the physics, 
engineering and mathematical literature. 

In the mathematical literature, CLTs for Wigner matrices can be traced back to Girko 
[9] (see also Q2]). Results for this class of matrices have also been obtained by Khorunzhy 
et al. [21j . Johansson [TH], Sinai and Sochnikov [25], Soshnikov [53], Cabanal-Duvillard [7]. 
For band matrices, let us mention the paper by Khorunzhy et al. |21j , Boutet de Monvcl 
and Khorunzhy [5], Guionnet [13], Anderson and Zeitouni pQ. The case of Gram matrices 
has been studied in Jonsson [20] and Bai and Silverstein [3]. For a more detailed overview, 
the reader is referred to the introduction in [Tj. In the physics literature, so-called replica 
methods as well as saddle-point methods have long been a popular tool to compute the 
moments of the limiting distributions related to the fluctuations of the statistics of the 
eigenvalues. 

Previous results and methods have recently been exploited in the engineering literature, 
with the growing interest in random matrix models for wireless communications (see the 
seminal paper by Telatar [30] and the subsequent papers of Tse and co-workers [31], [32] ; 
see also the monograph by Tulino and Verdu [33] and the references therein). One main 
interest lies in the study of the convergence and the fluctuations of the mutual information 
jj logdet (Y n Y* + pIn) for various models of matrices Y n . General convergence results have 
been established by the authors in [I7J [TS] [TB] while fluctuation results based on Bai and 
Silverstein [3] have been developed in Debbah and Muller [8] and Tulino and Verdu [34] . 
Other fluctuation results either based on the replica method or on saddle-point analysis 
have been developed by Moustakas, Sengupta and coauthors [22] [23], Taricco [29]. In a 
different fashion and extensively based on the Gaussianity of the entries, a CLT has been 
proved in Hachem et al. |14] . 

Comparison with existing work. There are many overlaps between this work and other works 
in the literature, in particular with the paper by Bai and Silverstein [3] and the paper by 
Anderson and Zeitouni [I] (although this last paper is primarily devoted to band matrix 
models, i.e. symmetric matrices with a symmetric variance profile). The computation of the 
variance and the obtention of a closed-form formula significantly extend the results obtained 
in [H]. 

In this paper, we deal with complex variables which are more relevant for wireless com- 
munication applications. The case of real random variables would have led to very similar 
computation, the cumulant k = E|X| 4 — 2 being replaced by k = EX 4 — 3. Due to the 
complex nature of the variables, the CLT in [T] does not apply directly. Moreover, we sub- 
stantially relax the moment assumptions concerning the entries with respect to [1] where the 
existence of moments of all order is required. In fact, we shall only assume the finitcness of 

the 8 moment. Bai and Silverstein [3] consider the model T„ 2 X n X*Tn where the entries 
of X n are i.i.d. and have gaussian fourth moment. This assumption can be skipped in our 
framework, where a good understanding of the behaviour of the diagonal individual entries 
of the resolvent (~zl n + YnY*)" 1 enables us to deal with non-gaussian entries. 

On the other hand, it must be noticed that we establish the CLT for the single functional 
logdet(F n y^ +pIn) and do not provide results for a general class of functionals as in [1] and 
[3]. We do believe however that all the computations performed in this article are a good 
starting point to address this issue. 
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Outline of the article. 

Non- asymptotic vs asymptotic results. As one may check in Theorems 13.11 13.21 and 13.31 
we have deliberately chosen to provide non-asymptotic (i.e. depending on n) determinis- 
tic formulas for the variance and the bias that appear in the fluctuations of X„(p). This 
approach has at least two virtues: Non-asymptotic formulas exist for very general vari- 
ance profiles (oij(n)) and provide a natural discretization which can easily be implemented. 
In the case where the variance profile is the sampling of some continuous function , i.e. 
o~ij(n) = a(i/N,j/n) (we shall refer to this as the existence of a limiting variance profile), 
the deterministic formulas converge as n goes to infinity (see Section [4]) and one has to 
consider Frcdholm determinants in order to express the results. 

The general approach. The approach developed in this article is conceptually simple. The 
quantity I n (/?) — EI„(p) is decomposed into a sum of martingale differences; we then system- 
atically approximate random quantities such as quadratic forms x T Ax where x is some ran- 
dom vector and A is some deterministic matrix, by their deterministic counterparts ^ Trace A 
(in the case where the entries of x are i.i.d. with variance i) as the size of the vectors and the 
matrices goes to infinity. A careful study of the deterministic quantities that arise, mainly 
based on (deterministic) matrix analysis is carried out and yields the closed-form variance 
formula. The martingale method which is used to establish the fluctuations of I n (p) can 
be traced back to Girko's REFORM (REsolvent, FORmula and Martingale) method (see 
El 112) ancl is close to the one developed in [3J. 

Contents. In Section^ we introduce the main notations, we provide the main assumptions 
and we recall all the first order results (deterministic approximation of WE n (p)) needed in 
the expression of the CLT. In Section [31 we state the main results of the paper: Definition of 
the variance 8^ (Theorem 13. 1[) , asymptotic behaviour (fluctuations) of N (T n (p) — EI?i(/o)) 
(Theorem I3.2[) . asymptotic behaviour (bias) of N (ET„(p) — Jlog(p + t)ir n (dt)) (Theorem 
13.3ft . Section [5] is devoted to the proof of Theorem 13. 1\ Section to the proof of Theorem 
13.21 and Section [3 to the proof of Theorem 13.31 

Acknowlegment. This work was partially supported by the Fonds National de la Science 
(France) via the ACI program "Nouvelles Interfaces des Mathematiques" , project MALCOM 
n° 205. 

2. Notations, assumptions and first order results 

2.1. Notations and assumptions. Let N = N(n) be a sequence of integers such that 

N(n) 

lim — = c e (0, oo) . 

n — >oo ri 

In the sequel, we shall consider a N X n random matrix Y n with individual entries: 

yn = g u( n ) X n 
ij ^ ij > 

where A™- are complex centered i.i.d random variables with unit variance and (cry(n); 1 < 
i < jV, 1 < j < n) is a triangular array of real numbers. Denote by var(Z) the variance of 
the random variable Z. Since var(Y^") = o~fj(ri)/n, the family (er^ (n)) will be referred to as 
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a variance profile. 



The main assumptions. 

Assumption A-l. The random variables (JT™ ; l<i<N,l<j<n,n>l) are complex, 
independent and identically distributed. They satisfy 

EX% = EpC™) 2 = 0, E|X™| 2 = 1 and E\X%\ 8 < oo . 

Assumption A-2. There exists a finite positive real number cr max such that the family of 
real numbers (<7jj(n), 1 < i < N, 1 < j ' < n, n > 1) satisfies: 

sup max \<Tij{n)\ < cr max . 

n>l )i% N 

— l<j<n 

Assumption A-3. There exists a real number cr 2 ^ > such that 

N 



lim inf min — cr 2 ,- in) > a 

n>l Ki<n n ^ 3 



2 

mill 



Sometimes we shall assume that the variance profile is obtained by sampling a function 
on the unit square of R 2 . This helps to get limiting expressions and limiting behaviours (cf. 
Theorem [231): 



Assumption A-4. There exists a continuous function a 2 : [0, 1] x [0, 1] — ► (0, oo) such that 
r? 3 ( 



aUn) = a 2 (i/N,j/n). 



Remarks related to the assumptions. 

(1) Using truncation arguments a la Bai and Silverstein [2j [24} [25] , one may lower the 
moment assumption related to the X+j's in A-fTI 

(2) Obviously, assumption A-[3]holds if afj is uniformly lower bounded by some nonneg- 
ative quantity. 

(3) Obviously, assumption A-|4] implies both A-[2]and A-[3] When A-[4] holds, we shall 
say that there exists a limiting variance profile. 

(4) If necessary, assumption A-[3]can be slightly improved by stating: 

( 1 N 1 n , \ 

max lim inf min — > o 2 An) , lim inf min — > cnAn) > . 

I n>l l<j<ri n .t-^ 3 n>\ Ki<N n ^ * 3 v ' I 

In the case where the first liminf is zero, one may notice that logdct(Y"„y^ + pl/v) = 
logdet(F„*y„ + pl n ) + (n — N) \ogp and consider Y*Y n instead. 

Notations. The indicator function of the set A will be denoted by 1a(x), its cardinality by 
#A As usual, l+ = {i£l : x > 0} and C+ = {z € C : ha(z) > 0}. 

V ... . T> 

Wc denote by — > the convergence in probability of random variables and by — > the 
convergence in distribution of probability measures. 

Denote by diag(ai; 1 < i < k) the k x k diagonal matrix whose diagonal entries are 
the ttj's. Element of matrix M will be either denoted rriij or [Af]y depending on the 
notational context. Denote by M T the matrix transpose of M, by M* its Hcrmitian adjoint, 
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by Tr(M) its trace and det(M) its determinant (if M is square), and by F M M , the empirical 
distribution function of the eigenvalues of MM*, i.e. 

F MM \x) = ^#{i:\ i <x} , 

where M M* has dimensions N x N and the Ai's are the eigenvalues of MM*. 

When dealing with vectors, | • | will refer to the Euclidean norm, and || • ||oo, to the 
max (or H.^) norm. In the case of matrices, || • || will refer to the spectral norm and 
III 'loo ^° ^ ne maximum row sum norm (referred to as the max-row norm), i.e., |||M||| = 
maxi<i<7v X^Li |[-^]u l when M is a N x N matrix. We shall denote by r(M) the spectral 
radius of matrix M . 

When no confusion can occur, we shall often drop subscripts and superscripts n for read- 
ability. 

2.2. Stieltjes Transforms and Resolvents. In this paper, Stieltjes transforms of proba- 
bility measures play a fundamental role. Let v be a bounded non-negative measure over R. 
Its Stieltjes transform / is defined as: 

/(*)= / *€C\supp(i/) , 

where supp(^) is the support of the measure v. We shall denote by <S(R + ) the set of Stieltjes 
transforms of probability measures with support in R + . 

We list in the following proposition the main properties of the Stieltjes transforms that 
will be needed in the paper: 

Proposition 2.1. The following properties hold true. 

(1) Let f be the Stieltjes transform of a probability measure v on R ; then: 

- The function f is analytic over C \ supp(^). 

- If f(z) G S(R + ), then \ f(z)\ < (d(z,R + ))" 1 where d(z,R+) denotes the dis- 
tance from z to R + . 

(2) Let P„ and P be probability measures over R and denote by f n and f their Stieltjes 
transforms. Then 

(VzeC+, f n (z) >f(z)) P„-^-,P. 

There are very close ties between the Stieltjes transform of the empirical distribution of 
the eigenvalues of a matrix and the resolvent of this matrix. Let M be a N x n matrix. The 
resolvent of MM* is defined as: 

Q{z) = (MM* - zIn)" 1 = (quiz))^^ , z e C - R+ . 

The following properties are straightforward. 

Proposition 2.2. Let Q(z) be the resolvent of MM*, then: 

(1) The function h n (z) = ^Tr Q(z) is the Stieltjes transform of the empirical distri- 
bution of the eigenvalues of MM* . Since the eigenvalues of this matrix are non- 
negative, h n (z) € 5(R + ). 
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(2) For every z £ C-R+, ||Q(z)|| < (d(z,R+)) _1 . In particular, ifp>0, \\Q(-p)\\ < 



2.3. First Order Results: A primer. Recall that T n {p) = -k log det (Y n Y* + pi) and let 
p > 0. We remind below some results related to the asymptotic behaviour of EJ n (p). As 

= M E l0 § ( A « +/')=/ l0 8'( A + ^ ( A ) , 

7V i=l 

where the Ai 's are the eigenvalues of YY* , the approximation of E2"„ (/?) is closely related to 
the "first order" approximation of F YnY ™ as n — > oo and iV/n — > c > 0. 

The following theorem summarizes the first order results needed in the sequel. It is a 
direct consequence of jTTJ Sections 2 and 4] (see also [TTjV 

Theorem 2.3 (|17j. Consider the family of random matrices (Y n Y*) and assume that 

A-Q] and A-[2] hold. Then, the following hold true: 

(1) The system of N functional equations: 

U(z) = to (2.1) 

z "t- „ Z^=i l+i^^w^w 

admits a unique solution (ti(z),--- ,t^(z)) in S(Bi + ) N . In particular, m n (z) = 
-k J^iLi belongs to S(M. + ) and there exists a probability measure ir n on R + such 
that: 

I \ f°° KnidX) 

m n (z) — 



X-z 

(2) For every continuous and bounded function g on R + ; 

g{\) dF Y " Y - (X) - ( g(X)7r n (d\) >0 a.e. 

R+ Jr+ n ^°° 



(3) The function V n (p) = L + log(A + p)ir n (dX) is finite for every p > and 
El n {p) - V n {p) ► where l n (p) = i- log det [Y n Y* + P I N ) . 

n — ►oo Jv 

Moreover, V n (p) admits the following closed form formula: 
1 - 1 n f 1 N \ 

V ^p) = -]vE lo s^(-^ + ^E lo s 1 + -E4»(-p) 

i=l 3=1 \ f=l / 

1 \ - a-^(n)ti(-p) 

Nn 2—1 



l= l:N.j = l:n 1 + £ E«=l ^ (»*)*<(-p) 

where the ti 's are defined above. 



Theorem 12.31 partly follows from the following lemma which will be often invoked later 
on and whose statement emphasizes the symmetry between the study of Y n Y* and Y*Y n . 
Denote by Q n {z) and Q n (z) the resolvents of Y n Y* and Y*Y n , i.e. 

Q n (z) = (Y n Y r : - zIn)- 1 = Mz))^^ , 2GC-K+ 



p + 



s 
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Lemma 2.4. Consider the family of random matrices (Y n Y*) and assume that A-fTI and 
A-[2] hold. Consider the following system of N + n equations: 

( U Jz) = -, ; - ; \ for 1 < i < N 

I tj n (z) = —, : — — —y for 1 < j < n 

{ J ' V ' z(l+ilYD 3l „T„(z)) J 

where 

T n {z) = diag(t J ,„(z), 1 < i < N), f n {z) = diagfo- n (z), 1 < j < n) , 

Dj. n = diag(o-?- (n), 1 < i < N), A> = diag(<4 (n),l<j<n). 

Then the following holds true: 

(a) |17[ Theorem 2.4] This system admits a unique solution 

{tl,ni ' ' ' ? tN,ni tl,n, ' ' ' tn,n) € S(R + ) N+n 

(Jo) \17[ Lemmas 6.1 and 6.6] For every sequence U n of N x N diagonal matrices and 
every sequence U„ of n x n diagonal matrices such as sup n max U n \\ , || U n \\) < oo , 
the following limits hold true almost surely: 

lim Ilk ([/„ (Q n (z) - T n (z))) = VzeC-R+, 

n—*oo,N/n — >c IV 

lim -Tr(U, (q„(z) - f n (z))) = VzeC-l+. 

n-ioo,JV/n-ic n V V // 

In the case where there exists a limiting variance profile, the results can be expressed in 
the following manner: 

Theorem 2.5 ([B], |10j . |16| ). Consider the family of random matrices (Y n Y*) and assume 
that A-Q] and A-|4] hold. Then: 

(1) T/ie functional equation 

*•«>- f-* + r 14 . .y, w *r < 2 - 2 » 

y Jo 1 + c J c (x, v)t(x, z) ax y 

admits a unique solution among the class of functions $ : [0, 1] x C \ R — > C smc/i 
£/ia£ u i— > z) is continuous over [0, 1] and z i— > <E>(u, z) belongs to S(JR + ). 

(2) T/ie function f(z) = J Q t(u, z) aw where t(u, z) is defined above is the Stieltjes trans- 
form of a probability measure P. Moreover, we have 

F r n r: P a . s . 



Remark 2.1. If one is interested in the Stieltjes function related to the limit of F Y ™ Yn , then 
one must introduce the following function f . which is the counterpart of r: 

f(v,z) = (— z + c [ ; 1-!—} dt.\ 

V J° l + fo<T 2 (t,s)T(s,z)ds J 
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Functions r and f are related via the following equations: 

t(u, z) = — ^- and f(v, z) = — - ^- . 

z I 1 + J <j 2 (u, v)t(v, z) dv j zil + cj a 2 (t, v)r(t, z) dt J 

(2.3) 

Remark 2.2. We briefly indicate here how Theorems 12.31 and 12.51 above can be deduced from 
Lemma l2~4l As ^-TrQ„(z) is the Stieltjes transform of F YnY " , Theoreml2~4l-f6) with U n = In 
yields -^-TrQ n (z) — -^-TrX^z) — + almost surely. When a limit variance profile exists as 
described by A-|4l one can easily show that j^TiT n (z) converges to the Stieltjes transform 
f(z) given by Theorem 12.51 (Equation (|2 . 2[) is the "continuous equivalent" of Equations 
(|2.ip ). Thanks to Proposition 12 . II - (|2|) . we then obtain the almost sure weak convergence of 
pY n Y n -j- j n .j^g cage w j iere A-|4]is not satisfied, one can prove similarly that F YnY ™ is 
approximated by 7r„ as stated in Theorem 12 . 31 - ( 2) . 



3. The Central Limit Theorem for I n {p) 

When given a variance profile, one can consider the V s defined in Theorem l2.3K l). Recall 
that 

T(z) = diag(*i(z),l < i < N) and D } ■= diag(cr J 2 J „ 1 < i < N) . 

We shall first define in Theorem 13.11 a non-negative real number that will play the role of 
the variance in the CLT. We then state the CLT in Theorem 13.21 Theorem 13.31 deals with 
the bias term N(KZ - V). 

Theorem 3.1 (Definition of the variance). Consider a variance profile (o~ij) which fulfills 
assumptions A-121 and A-[3] and the related ti 's defined in Theorem \2.3\ -fl). Let p > 0. 



(1) Let A n = [ai^m) he the matrix defined by: 

1 ±-TrD e D m T(-p) 2 



1 < I , m < n 



n(l + lTvD e T(-p)y 
then the quantity V n = — logdet(/ n — A n ) is well-defined. 

(2) Denote by W n = Tryl n and let k be a real numbed satisfying k > — 1. The sequence 
(V n + fcW n ) satisfies 

< liminf (V„ + kW„) < limsup (V„ + kW„) < oo 

" n 

as n — > oo and N/n — > c > 0. We shall denote by: 

01 = - logdet(J - A n ) + kTtA u . 



Proof of Theorem 13.11 is postponed to Section [5j 

In the sequel and for obvious reasons, we shall refer to matrix A n as the variance matrix. 
In order to study the CLT for N(T n (p) — V n {p)), we decompose it into a random term from 
which the fluctuations arise: 

n (i n { P ) - m n (p)) = io g det(r n y n * + P i N ) - Eiogdet(y n y„* + P i N ) , 



2 In the sequel, k is denned as k = E|Xn| 4 — 2. 
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and into a deterministic one which yields to a bias in the CLT: 

N(¥J n (p) - V n {p)) =mogdct{Y n Y:+pI N ) -N J log(A + p)% n (d\) . 
We can now state the CLT. 

Theorem 3.2 (The CLT). Consider the family of random matrices (Y n Y*) and assume 
that A-HJ A-[2] and A-[3] hold true. Let p > 0, let n = E|Xn| 4 - 2, and let 0^ be given by 
Theorem \3.1l Then 

e,- 1 ( log det(y„y„* + P i N ) - e log dct(r„y;* + P i N )) - > M(o, 1) . 



Proof of Theorem 13.21 is postponed to Section [5] 

Remark 3.1. In the case where the entries are complex Gaussian (i.e. with independent 
normal real and imaginary parts, each of them centered with variance 2 _1 ) then k = and 
2 reduces to the term V n . 



The asymptotic bias is described in the following theorem: 
Theorem 3.3 (The bias). Assume that the setting of Theorem VS .'A holds true. Then 

(1) For every lj G [p,+oo), the system of n linear equations with unknown parameters 
(w e . n (uj); l<£<n): 

m=l ^ n * V /' 



H = ^ 2 t £ (-c) 2 (^y f4Mz^! T rA 2 T(-c) 2 ) - M^) Tr ^ 2 T(-^ 2 



n 

i=l 



(3.2) 

admits a unique solution for n large enough. In particular if k = 0, then p e n = 
and to^ iTt = 0. 
(2) Let 

1 ™ 

/3„(w) = - V"w^„(w) . (3.3) 
n ' 

Then B n (p) = J°° f3 n (u>) duj is well-defined, moreover, 

poo 

limsup / \(3 n {u))\du) < oo. (3-4) 
n J p 

Furthermore, 

N(El n (p)-V n (p))-B n (p) ■ ► 0. (3.5) 



Proof of Theorem 13.31 is postponed to Section [7] 
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4. The CLT for a limiting variance profile 

In this section, we shall assume that Assumption A-|4] holds, i.e. cr^-(n) = <r 2 (i/N, j/n) 
for some continuous nonnegative function a 2 (x,y). Recall the definitions (|2.2p of function 
r and of the tj's (defined in Theorem 12.31 (1)). In the sequel, we take p > 0, z = —p and 

denote r(t) = r(t, —p). We first gather convergence results relating the t^s and r. 

Lemma 4.1. Consider a variance profile (cry) which fulfills assumption A-|4j Recall the 
definitions of the ti 's and r. Let p > and let z = —p be fixed. Then, the following 
convergences hold true: 

(1) W ti& i — - — > T ( u ) du , where ^> stands for the weak convergence of measures. 

N n— >oo 

(2) sup i<N \ti-T(i/N)\ ,0 . 

— n — 'oo 

(3) i^nS^^rH^du, 

Proof. The first item of the lemma follows from Lemma I2.4l fb) together with Theorem 
2.3-(3) in [16]. 

In order to prove item (2), one has to compute 
1 " a 2 (i/N,j/n) 



U - r(i/N) 



i=l 1 + n 




a 2 (u,v) , > 

-dv 



la 1 + c (t 2 (x,v)t(x) dx 

and use the convergence proved in the first part of the lemma. In order to prove the unifor- 
mity over i < N, one may recall that C[0, l] 2 = C[0, 1] <g> C[0, 1] which in particular implies 
that Ve > 0, there exist gi and hg such that sup^. y \a 2 (x,y) — X^=i 9e( x )hi(y)\ — e - Details 
are left to the reader. 

The convergence stated in item (3) is a direct consequence of item (2). □ 

4.1. A continuous kernel and its Predholm determinant. Let K : [0, l] 2 — > R be 

some non-negative continuous function we shall refer to as a kernel. Consider the associated 
operator (similarly denoted with a slight abuse of notations) : 

K:C[0,1] -» C[0,1] 



/ » Kf(x)= / K(x,y)f{y)dy . 

J [0,1] 

Then one can define (see for instance [37J Theorem 5.3.1]) the Frcdholm determinant det(l + 
XK), where 1 : / i— > / is the identity operator, as 

det(l - A/v) = V tll^l [ K ( Xl Xk ) ®U dXi (4.1) 

where 



: x y \ 2 , 
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for every A € C. One can define the trace of the iterated kernel as: 



'[04] 

In the sequel, we shall focus on the following kernel: 

K 00 (x,y) = —^ LJ -2-. (4.2) 

( 1 + c J, Q - cr 2 (u, x)t(u) du) 

Theorem 4.2 (The variance). Assume that assumptions A-Q] and A-^hold. Let p > and 

recall the definition of matrix A n : 

l „2 

m < n 



1 iiZjH" VjV'rJ" UVj_nJ_! 1 < £ 

n (i + ^Ef =1 - 2 (i4)^ 



Then: 

(1) TrA„ ► Tr^ . 

n — >oc 

(2) det(/„ - A n ) ► dct(l - K^) and dct(l - K^) ^ 0. 

(3) Let « = E|X n | 4 - 2, i/ien 

< - logdet(l - A'oo) + kTtKoo < oo . 

Proof. The convergence of TrA„ toward Tri^oo follows from Lemma T4. 11 (1) .(3). Details of 
the proof are left to the reader. 

Let us introduce the following kernel: 

One may notice in particular that at_ rn = ^-K n (—, ^). Denote by || • the supremum norm 
for a function over [0, l] 2 and by cr^ax = 1 1 ' 2 1 loo, then: 

||*»||oo < and HJiCoolloo < c^f . (4.3) 

n p^ p z 

The following facts (whose proof is omitted) can be established: 

(1) The family (if„)„>i is uniformly equicontinuous, 

(2) For every (a;, y), K n (x, y) -> K^x, y) as n -> oo. 

In particular, Ascoli's theorem implies the uniform convergence of K n toward K^. It is now 
a matter of routine to extend these results and to get the following convergence: 

K n {^ Xk ) > Koo (^ Xk ) (4.4) 

V yi ■■■ Vk J «^oo V vi ■ • • Vk J 
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uniformly over [0, l] 2fc . Using the uniform convergence (|4.4|) and a dominated convergence 
argument, we obtain: 

- k £ kJ */» »/» ) — ► / K J* Xk ) ®u d*< ■ 

n * t-j V i x n ■■■ i k /n J n->oo J, 01]k \ x\ ••• x k I 

l<zi ,2fc<n 



Now, writing the determinant det(/„ + \A n ) explicitely and expanding it as a polynomial in 
A, we obtain: 

det(/„-AA n ) = £H^(% £ K n {^ !f 

fc=0 \ l<ii,— ,tj,<n v 

Applying Hadamard's inequality ( [271 Theorem 5.2.1]) to the determinants and -Koo( - ) 

yields: 



l<ii ,?fc <n 

where (a) follows from (|4.3[) . Similarly, 



^ , V *i/ n • • • l k/n J 



[0,1] 



K x [ Xl "' Xfc da* < 

1 xi • ■ • x fc 



Since the series ,f |A| converges, a dominated convergence argument yields the con- 

vergence 

det(7„ + AA„) ► det(l + , 

n — >oo 

and item (2) of the theorem is proved. Item (3) follows from Theorem 13. 11 ( 2) and the proof 
of the theorem is completed. 

□ 



4.2. The CLT: Fluctuations and bias. 

Corollary 4.3 (Fluctuations). Assume that (A\J\) and (A^ hold. Denote by 

9^ = - logdet(l - A'oo) + KTrK^ , 

then 

N 

— (X n (p)-EX„(p)) 

= G^ 1 (log det {Y n Y* + P I N )-E log det {Y n Y* + P I N )) — ±-> Af(0, 1) . 

n — >oo 

Proof, follows easily from Theorem 13.21 and Theorem 14.21 □ 
Recall the definition of f (cf. Remark |2.1[) . 
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Theorem 4.4 (The bias). Assume that the setting of Corollary \4-3\ holds true. Let oj G 
[p, oo) and denote by p : [0, 1] — > R the quantity: 

p{x,U)) — KOJ 2 T 2 (x,—Ld) 

x jwc J a 2 (u,x)T Z {u) (^J a 2 (s,u)f 2 (s)ds^j du 

T/ie following functional equation admits a unique solution: 

i \ /" 1 c/ 1 (T a (u,a)tr a (u,y)r a (u)du u , , > 

w{x,uj)= \ -5-iu(y,w)dj/ + p(x,w) . 



/o (l + C r o 1 C r2( U;2 .) T ( u ) du y 

Lei ^ooC^) = Jq w(x,u>) dx . Then J°° |/?oo(^)| < °°- Moreover, 

roc 

N{m n {p)-V n {p)) ► B OQ {p)= /?ooM^. (4.5) 



Proof of Theorem 14.41 although technical, follows closely the classical Fredholm theory as 
presented for instance in [2TJ Chapter 5]. We sketch it below. 

Sketch of proof. The existence and unicity of the functional equation follows from the fact 
that the Fredholm determinant det(l — Koo) differs from zero. In order to prove the con- 
vergence (|4.5[) . one may prove the convergence J f3„ — > J /3oo (where 0„ is defined in 
Theorem 13. 3p by using an explicit representation for /3oo relying on the explicit representa- 
tion of the solution w via the resolvent kernel associated to K x (see for instance [27l Section 
5.4]) and then approximate the resolvent kernel as done in the proof of Theorem 14.21 □ 

4.3. The case of a separable variance profile. We now state a consequence of Corollary 
14.31 in the case where the variance profile is separable. Recall the definitions of r and f given 
in (TO). 



Corollary 4.5 (Separable variance profile). Assume that A-[T] and A-l4l hold. Assume more- 
over that p > and that a 2 is separable, i.e. that 

a 2 (x,y) = d(x)d(y) , 

where both d : [0, 1] — ► (0, oo) and d : [0, 1] — > (0, oo) are continuous functions. Denote by 



l j-i 



7 = c/ d 2 (t)T Z (t)dt and 7= / d z (t)f 2 (t) dt . 
Jo Jo 

Then 

8^ = -log(l-p 2 7 7 ) +^77 . (4.6) 

Remark 4.1. In the case where the random variables are standard complex circular 
gaussian (i.e. Xij = Uij + iVij with [Zy and V^j independent real centered gaussian random 
variables with variance 2 _1 ) and where the variance profile is separable, then 



TV (0, - log ( 1 

— *oo 

This result is in accordance with those in [22] and in [T^j . 



iV(X„(p) - V n (p)) — ^ A/" (0, - log (1 - P 2 77)) 

n — *-oo 
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Proof. Recall the definitions of r and f given in (|2.3[) . In the case where the variance profile 
is separable, the kernel writes: 

cd(x)d(y) J, Q j, d 2 {u)T 2 (u) du _ 
Koo(x,y) = ! -2 = p jd(x)d(y)f (x) . 



1 + cd(x) Jj Q ^ d(u)r(u) du 

In particular, one can readily prove that TrA'oo = p 2 ^. Since the kernel K QO (x 1 y) is itself 
a product of a function depending on x times a function depending on y, the determinant 

Kryo (^ x y \ ... ) i s equal to zero for k > 2 and the Frcdholm determinant writes det(l — 
Koo) = 1 - Jjo.i] Koc(x,x)dx = 1 - p 2 77 . This yields 

- logdet(l - A 00 ) + kTtKoo = - log(l - ^77) + Kp 2 ^ > 
which ends the proof. □ 

5. Proof of Theorem 13. II 

Recall the definition of the n x n variance matrix A n : 

1 TvD e D m T(-p) 2 
n 2 (1 + ±TrD e T(-p)) 

In the course of the proof of the CLT (Theorem 13. 2j) , the quantity that will naturally pop 
up as a variance will turn out to be: 

e 2 n = v„ + «w„ (5.1) 

(recall that W n = TrA„) where V n is introduced in the following lemma: 

Lemma 5.1. Consider a variance profile (cry) which fulfills assumptions A-121 and A-[3] and 
t/ie related ti 's defined in Theorem \2.!A (l). Let p > and consider the matrix A n as defined 
above. 

(1) For 1 < j < n, the system of (n — j + 1) linear equations with unknown parameters 



{y { H 3<e<n): 



Vl,l= E °<.m V#!n + (5-2) 
admits a unique solution for n large enough. 
Denote by V„ the sum of the first components of vectors (y^ j < t < n), i.e.: 

3=1 

2. Let k be a real number satisfying k > — 1. 27ie sequence [v n + KW m j satisfies 
< liminf ( V„ + kW„ ) < limsup ( V„ + kW„ ) < 00 



osmco and iV/n — > c > 0. 
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3. The following holds true: 

V n + logdet(I„ - A n ) 



0. 



Obviously, Theorem 13.11 is a by-product of Lemma 15.11 The remainder of the section is 
devoted to the proof of this lemma. 

We cast the linear system (|5.2[) into a matrix framework and we denote by A\P the 

(n - j + 1) X (n - j + 1) submatrix = (a e , m )7. m=J , by the (n - j + 1) x (n - j + 1) 

matrix A^) where the first column is replaced by zeros. Denote by the (n — j + 1) X 1 
vector : 

1 iTr^T(-p) 2 \ 
n 



(i+iw(-p)r;^ 

These notations being introduced, the system can be rewritten as: 

y£ ) =< ' ) 3/£' ) +<#' ) ^ (J->#tf>)i#>=d£>. (5.3) 

The key issue that appears is to study the invertibility of matrix (I — An ) and to get some 
bounds on its inverse. 



5.1. Results related to matrices with nonnegative entries. The purpose of the next 
lemma is to state some of the properties of matrices with non-negative entriess that will 
appear to be satisfied by matrices An . We shall use the following notations. Assume that 
M is a real matrix. By M y (resp. M y 0) we mean that my > (resp. my > 0) for 
every element my of M. We shall write M y M' (resp. M y M') if M - M' y (resp. 
AI — AI' y 0). If x and y are vectors, we denote similarly x y 0, x y and x y y. 

Lemma 5.2. Let A = (af,m)™ jra=1 &e a n x n real matrix and u — (u£, 1 < £ < n), 
v = (v£, 1 < £ < n) be two real n x 1 vectors. Assume that A y 0, u y ; and v y 0. 
Assume furthermore that equation 



u 



Au + v 



is satisfied. Then: 

(1) The spectral radius r(A) of A satisfies r(A) < 1 — < 1- 

(2) Matrix I n ~ A is invertible and its inverse (I n — A) 1 satisfies: 



(In 



A) 1 y0 and 



(In - Ay 



> 1 



/or every 1 < £ < n. 

(3) T7ie max-row norm of the inverse is bounded: (I n — A) 1 

(4) Consider the (n — j + 1) X (n — j + 1) submatrix A^' = (a£ m ) r t 



< 



min^(^) 



f and denote by 

A0,U) matrix A^' whenever the first column is replaced by zeros. Then properties 
(1) and (2) are valid for A® and 



I 



(n—j + l) 



< 



maxi<K ffi (Mf) 
mini<£<„(^) 



Proof. Let a = 1 



mm ^, . Since u y and v y 0, a readily satisfies a < 1 and ait ^= u — v = 
Au which in turn implies that r(A) < a < 1 jT5J Corollary 8.1.29] and (1) is proved. In order 
to prove (2), first note that Vra > 1, A m y 0. As r(A) < 1, the series J2 m >o^ m converges, 
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matrix /„ — A is invertible and (/„ — A) 1 = J2 m >o A m In fc= 0. This in particular implies 
that [(/„ — ^4) _1 ]« > 1 for every 1 < £ < n and (2) is proved. Now u = (I n — A)~ 1 v implies 
that for every 1 < k < n, 



Uk 



n n 

[(In - Ay 1 ] ^ Vl > min(^) [(^ ~ AT 1 

e=i ke i=i 



hence (3). 

We shall first prove (4) for matrix A^\ then show how A '^ inherits from A^'s proper- 
ties. In [18] , matrix A^ is called a principal submatrix of A. In particular, r(A^) < r(A) 
by [HI Corollary 8.1.20]. As A )p 0, one readily has A^ )p which in turn implies prop- 
erty (2) for Aw. Let A^ be the matrix A»' augmented with zeros to reach the size of 
A. The inverse (/„_j_|_i — A^)~ x is a principal submatrix of (J„ — A^')~ l )p 0. There- 



fore, 



[I(n-j+l) 



< 



(in-A^y 



Since A m ^ (A^) m for every m, 



one has £ m>0 A m > Y, m >o(A U) ) m ; equivalently (I - A)- 1 (I - A^)' 1 which yields 



< 



(i -Ay 



Finally (4) is proved for matrix A^> 



We now prove (4) for A°^\ By pH Corollary 8.1.18], r(A '^) < r(A^) < 1 as A® )>= 
A°'0). Therefore, (I - A ^) is invertible and 

(I - A '^)- 1 =jt[A°' U) ] k . 

fe=0 

This in particular yields (I - A '^)' 1 )? and (I - A°>^)^ > 1. Finally, as A® > A°^\ 
one has 



Item (4) is proved and so is Lemma. 15.21 



< 



i - A® 



□ 



5.2. Proof of Lemma I5.lt Some preparation. The following bounds will be needed: 

Proposition 5.3. Let p > 0, consider a variance profile (o~ij) which fulfills assumption A-121 
and consider the related ti 's defined in Theorem \2.!A (l). The following holds true: 



> U(-p) > 



1 



Proof. Recall that t e (z) E 5(R+) by Theorem O In particular, t e (-p) = J R+ t^jffi for 
some probability measure p(. This yields the upper bound tt(—p) < p^ 1 and the fact that 
U(—p) > 0. Now the lower bound readily follows from Eq. (|2.ip . □ 

Proposition 5.4. Let p > 0. Consider a variance profile ((Jy) which fulfills assumptions 
A-[2] and A-[3l' consider the related U 's defined in Theorem \2.3\ -fl). Then: 

liminf min — TrZLTL (— p) 2 > and liminf min —TrD 2 T n (—p) 2 > . 

n>l Kj<n n ra>l Kj<ra 71 
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Proof. Applying Proposition 15. 31 yields: 

which is bounded away from zero by Assumption A-[3l Similarly, 

which remains bounded away from zero for the same reasons (notice that (a) follows from 
the elementary inequality (n" 1 ^2 x i) 2 < n" 1 ~^2 X 1)- D 



We are now in position to study matrix A n = A n . 

Proposition 5.5. Let p > 0. Consider a variance profile (cry) which fulfills assumptions 
A-[2] and A-[3l' consider the related ti 's defined in Theorem \2.3\ -fl ) and let A„ be the variance 
matrix. Then there exist two n x 1 real vectors u„ = (u£ n ) y and v n = (vg n ) >- such that 
u n = A n u n + v n . Moreover, 

sup max (uf n ) < oo and liminf min (v» n ) > . 

n K£<n n l<£<n 



Proof. Let z = — p + <5i with 8 G K — {0}. An equation involving matrix A n will show up by 
developing the expression of Im(T(z)) = (T(z) — T*(z)) /2i and by using the expression of 
the tj(z)'s given by Theorem 12. 3K 1). We first rewrite the system (|2.1| as: 



TrD m T i 

\ m— l ™ / 

We then have 

Im(T) = i_( T _ T *) = i_ TT * |^ r *-i _ r - 

1 A D m TT* fl 

rim -TrD m T + STT 



m=l | 1 ^ n i - ylJ m ± | 
This yields in particular, for any 1 < £ < n: 



\m ( —Ti D^t] = \ V ^^£E1 i lm ( I^t) + —Tr Z^TT* . (5.5) 
V" / n ^ \l + ±TTD m T\ s \ n J n 



S 



Recall that for every 1 < i < N, ti(z) G 5(K + ). Denote by pi the probability measure 
associated with U i.e. ti(z) — J R+ . Then 

\ lm (l TrDlT \ = ifv, f°° J^L — »iv 



Denote by U£„ the right handside of the previous limit and let ug n = ^ 1+ i Tr ^" r ^_ p ^ 2 
Plugging this expression into (|5.5p and letting <5 — > 0, we end up with equation: 

u n A n u n -\- v n 
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A n )p is given in the statement of the lemma, u n = (u^ n ; 1 < £ < n) and v n = (v£ n ; 1 < 
I < n) are the n x 1 vectors with elements 



ue.; 



(1 + ±TrD e T(-p)) 



— — and V£ t „ — - 



{l + ±TrD e T(-p)y 



(5.6) 



It remains to notice that u n >- and u n >~ for n large enough due to A-[3j that the 
numerator of ut >n is lower than (7Vcr^ lax )/(np 2 ) and that its denominator is bounded away 
from zero (uniformly in n) by Propositions 15.31 and 15.41 Similar arguments hold to get a 
uniform upper bound for V(, n . This concludes the proof of Proposition 15.51 □ 

5.3. Proof of Lemma I5.lt End of proof. 

Proof of Lemma [5J\ -(l ). Proposition 15 . 51 together with Lemma [572]- (4) yield that I — A°'V' 
is invertible, therefore the system [5.31 admits a unique solution given by: 



y 



(i) = (/-4^))-id<?), 



and (1) is proved. 



□ 



Proof of Lemma \5.1Y (2). Let us first prove the upper bound. Proposition 15.51 together with 
Lemma 15.21 yield 



limsupmax (/ — A '^) 1 



maxi<£< n (i^„) 
< lim sup — : — - < oo 

„>i mmi<<<„(!)( n ) 



Each component of vector d^) satisfies e?£^ < N °2™i x i.e. sup 1<J<ri ||<i^||oo < lr- There- 



fore, vector yg) satisfies: 

BUpllt/^Hco < BUpKJ-^J-^oolld^lloo < £ ■ 
j j n 

Consequently, 

n n 

o < v„ = $>g ^ Eh^iu 

3=1 3=1 

satisfies lim sup n V„ < oo. Moreover, Proposition O yields VV„ < n~ 2 £? =1 Tr£> 2 T 2 < 
a^ lax N(p 2 n)~ 1 . In particular, W„ is also bounded and limsup n (V n + kW„) < lim sup ra ( V n + 
\n\W n ) < oo. 

We now prove the lower bound. 



Recall that = (/- A°'^)~ 1 d < £> . We therefore have: 



,,U) _ 



((/-^o.W))-i_/) d«) 



'(J-^W)- 1 ^^) 
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As (I - A°' (i) )- X > I, we have: 



Mi,™ 3,n — 



40,(j) d (j) 



E 



1 {^D t D 3 T\-p))' 



^ (l + ^ DjT {-p)) 



> * E - ( r Tr ^^ j 



where (a) follows from Proposition 15. 31 which is used both to get a lower bound for the 
numerator and an upper bound for the denominator: (1 + ^-TrDjT) 4 < (l+Nn~ 1 (T^ aax p~ l ) i . 
Some computations remain to be done in order to take advantage of A-[3] to get the lower 
bound. Recall that — X^fcLi x 1 — (m ^2T=i x k) ■ We have: 



E * EE 



3 = 1 l=j+l 

1 n(n - 1) 

n X " X 



n(n — 1) 



ttEO^ 



(a) 1 
> - 



(b) 1 

3 



> 



> 



> - 



by E ^Wi) +0 

l<j,f<n / 




0(1) 



«=1 J=l 



AT 



n ^ JV 

Ea7E^ 

,J = 1 i=l 




0(1) 



0(1) 



where (a) follows from the bound > | valid for n large enough. The term o(l) at 

step (6) goes to zero as n — * oo and takes into account the diagonal terms in the formula 
^ ^2j<£ a i e Tlj a iJ = t a j£- ^ rem a m s now to take the liminf to obtain: 



Item (2) is proved. 



□ 
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Proof of Lemma \5.1Y (3). We first introduce the following block-matrix notations: 



A® 



d91 <#> 



We can now express the following inverse: 



1 

(I — A 



and A2: U) 



{I-A^r 1 



a 



0) 



A 



(i+i) 



1 ^(/-A^)" 1 



(I — A 



O'+ih-i 



This in turn yields = dj^ + a 

vf\ < ~j where K does not depend on j and n, as 



(I — An ) 1 dffl and one can easily check that 



'3,n — n 



Remark that 



lyjSl^nl+^IU (i-A^r 1 ||d«)|| 



logdet(/ - A^) - logdet(J - 

=0') 



= log det 



log det 



1 - d}'^ -az 



(J-A^) 



(j'+lh-l 



1-d^ -^(/-^V 1 



-d# 



= log 1 - d 



(i) 



a^)(J-^+ 1 ))- 1 ^)) 



and write logdet(I — A n ) as: 

n-l 

logdet(/ -A n ) = E (logdct(7 - A&) - log det (J - + log(l - a nn ) 

i=i 

n-l 

= E log (l - dg» - - ^' +1) )- 1 di i) ) + log(l - a„„) 

n-l 

E (4S + &v - ^r 1 ^) + o (i) 

.to 



n-l 
-Vn + 0(1) • 



E^+°« 



This concludes the proof of Lemma 15.11 



□ 



6. Proof of Theorem 13.21 



6.1. More notations; outline of the proof; key lemmas. 
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More notations. Recall that Y„ — (y™) is a JV x n matrix where y™ = ^=X,- 7 - and that 
Qn(z) = {qij{z) ) = (Y n Y* - zI N )-\ We denote 

(1) by Q n (z) = (g y (*) ) = (Y*Y n - zI n )-\ 

(2) by yj the column number j of Y n , 

(3) by Y£ the iV x (n — 1) matrix that remains after deleting column number j from y n , 

(4) by Qj t n{z) (or Qj(z) for short when there is no confusion with Q n (z)) the N x N 
matrix 

Q J (z) = (y^Y^-z/ iV )- 1 , 

(5) by £i the row number i of Y„, 

(6) by Y^ n (or y for short when there is no confusion with Y n ) the (N — 1) x n matrix 
that remains after deleting row i from Y , 

(7) by Qi } n(z) (or Qi(z)) the n x n matrix 

= (r/y.-z/,,)- 1 . 

Recall that we use both notations qij or [<2],j for the individual element of Q{z) depending 
on the context (same for other matrices). The following formulas are well-known (see for 
instance Sections 0.7.3 and 0.7.4 in [TS]): 



" jj(z) = (6 - 2) 

For 1 < j < n, denote by Tj the cr-field Tj = cr(yj, ■ ■ ■ ,y n ) generated by the random 
vectors (%,••• ,y n )- Denote by Ej the conditional expectation with respect to !Fj, i.e. 
Ej = E(- | !Fj). By convention, J^ n +i is the trivial cr-field; in particular, E n +i = E. 

Outline of the proof. In order to prove the convergence of 0~ 1 (logdet(y„y jl * + plisr) — 
Elogdet(y„Y,* + pIn)) toward the standard gaussian law A/"(0, 1), we shall rely on the 
following CLT for martingales: 

Theorem 6.1 (CLT for martingales, Th. 35.12 in [3]). Let 7« , 7„_ 1 , . . . , 7i be a martin- 
gale difference sequence with respect to the increasing filtration Jn \ . . . , JF^ . Assume that 
there exists a sequence of real positive numbers 0„ such that 



T^E^r 2 -^!- (6-3) 
Assume further that the Lindeberg condition holds: 



n J=1 \ i 







TTien 0„ 1 X)j=i 7j"^ converges in distribution to A/"(0, 1). 
Remark 6.1. The following condition: 



2+5 
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known as Lyapounov's condition implies Lindcberg's condition and is easier to establish (see 
for instance [3], Section 27, page 362). 

The proof of the CLT will be carried out following three steps: 

(1) We first show that logdet(Y„F„* + pi) - Elogdet(y„F„* + pi) can be written as 

n 

io g det(y„y„* + P i) - Eio g det(y„y„* + P i) = 

3=1 

where (7^) is a martingale difference sequence. 

(2) We then prove that (7^) satisfies Lyapounov's condition (|6 .4|) where 2 is given by 
Theorem 13.11 

(3) We finally prove (|6 . 3[) which implies the CLT. 

Key Lemmas. The two lemmas stated below will be of constant use in the sequel. The first 
lemma describes the asymptotic behaviour of quadratic forms related to random matrices. 

Lemma 6.2. Let x = (x±, ■ ■ • ,x n ) be a n x 1 vector where the Xi are centered i.i.d. complex 
random variables with unit variance. Let M be a n x n deterministic complex matrix. 

(1) (Bai and Silverstein, Lemma 2.7 in j^jj Then, for any p > 2, there exists a constant 
K p for which 

E\x*Mx - TtM\ p < K p ((E|a;i| 4 TrA/M*) p/2 +E|.Ti| 2p Tr(MM*) p / 2N 



(2) (see also Eq. (1.15) in [5],) Assume moreover that Ex\ =0 and that M is real, the 

n 

E (x*Mx - TiM) 2 = TrM 2 + k ^ ; 
where k = E|xi| 4 — 2. 



i=l 



As a consequence of the first part of this lemma, there exists a constant K independent 
of j and n for which 

p 

< Kn- p l 2 (6.5) 



E 



VjQi(-p)vi-^DjQj(-p) 



for p < 4. 

We introduce here various intermediate quantities: 

d(z) = — — r , l<i<N; C(z) = diag( Cl (z); 1 < i < N) 

1 + iTrAEQ(z) 



&j(z) = — 7 3 — -T-, 1 < j < n: C(z) = diag(c,(z): 1 < 7 < n) , 

j\ ) z (l + ±Tr DjEQ(z)) - ■> - > ^ 51,1 -J- J> 

bi(z) = —. — : r , l<i<N; B(z) = diag(6 8 (z); 1 < i < N) 

1 + iTrAC(z) 



bj(z) = j- j-j- — v 1 < j < n; B(z) = diag(6,(z); 1 < j < n) . (6.6) 
z(l + ^TrDjCiz)) 

The following lemma provides various bounds and approximation results. 
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Lemma 6.3. Consider the family of random matrices (Y n Y*) and assume that A-fTI and 
A-[2] hold true. Let z = —p where p > 0. Then, 

(1) Matrices C n satisfy \\C n \\ < - and < Cj < - . These inequalities remain true 
when C is replaced with B or C . 

(2) Let U n and U n be two sequences of real diagonal deterministic N x N and nx n ma- 
trices. Assume that sup Jt>1 max M|f n ||, ||t^n||J < oo , then the following hold true: 

(a) -TrlI(EQ - T) > and -TrC7(EQ - f) > , 

(b) ±TrU(B - T) ► , 

n — >oo 

(c) sup„ E (TW (0 - EQjf < oo. 

(d) sup„ 4,E (TrU (Q - EQ)) 4 < oo, 

(3) [Rank-one perturbation inequality] The resolvent Qj satisfies \TrM (Q — Qj) \ < 



\M || 
P 



for any N x N matrix M (see Lemma 2.6 in [25].)- 
Proof of Lemma 16.31 is postponed to Appendix [5] 

Finally, we shall frequently use the following identities which arc obtained from the defi- 
nitions of c, and Cj together with Equations 



[Q{z)]u = Ci + zc4Q] u UiQi& - ^TtDiEQj (6.7) 
[Q(z)]jj = Cj+zcMn '(vjQjVj ■■-^DjEQj (6.8) 

6.2. Proof of Step 1: The sum of a martingale difference sequence. Recall that 
Ej = E(- | T 3 ) where T$ = a{y e , j <£< n). We have: 

io g dct(ry* + P i N ) - Eio g dct(yr* + P i N ) 

n 

= ^(Ej - e j+1 ) log dct(yr* + P i N ) 

(a) (det (YW + pI N )\ 



det (YY* + P I N ) 



= -2^( E J - e j+i) 1o s 



(c) 



(<*) 



n 

^(E, -E J+1 )log [Q(-p) 



det (F*F + p/„) 



2 (Ej - %+i) log (i + t£Qi (-p) %) 
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where (a) follows from the fact that does not depend upon yj, in particular Ej logdet {Y^Y' J + pl^) 
Ej+i logdet (YiY j * + pI N ) ; (b) follows from the fact that det(Y^*Y^ +pl n -i) = det{Y^Y j * + 
pIn) x p"" 1 "^ (and a similar expression for dct(Y*Y + pl n ))', (c) follows from the equality 

- _ dct(yj*yj + P j n _ 1 ) 
m pn "~ det(r*y + P /„) 

which is a consequence of the general inverse formula A^ 1 = dc ^ A ^ adj (^4) where adj(^l) is 
the transposed matrix of cofactors of A (see Section 0.8.2 in [IB]); and (d) follows from 
We therefore have 

io g dct(yy* + pi N )- Eio g dct(yy* + P i N ) 

n n 

= ^(Ej-E.+Olog^l + ^Q^-p)^) ^ £ 7j . 

As the following identity holds true, 

log M + -TrDjQA = E j+1 log M + -TrDjQj 



one can express 7j as: 

7j = (Ej - Ej+i) log [ 1 



1 + iTr^-Q, 



= (E J -E J+1 )log(l + r j ) where 1', ' V- -Tr/),(^ ' (6 ' 9) 

The sequence 7 „ , . . . , 71 is a martingale difference sequence with respect to the increasing 
filtration J- n , . .., T\ and Step 1 is established. 

6.3. Proof of Step 2: Validation of Lyapounov's condition (I6.4p . In the remainder 
of this section, z = —p. Let 5 > be a fixed positive number that will be specified below. 
As liminf 6£ > by Theorem [3TT1 we only need to prove that YJj=i E l7j \ 2+S 0- We 
have E|7j| 2+(5 = E \(Ej — E_j +1 ) log(l + r_j)| 2+<5 ; Minkowski and Jensen inequalities yield: 

(E| 7 ,| 2+i )^ < (E|E J log(l + r j )| 2+5 )^ + (E|E J+1 log(l + r j )| 2+5 )^ 
< 2(E|log(l + r j )| 2+5 ^ 



Otherwise stated, 

E| 7 ,| 2+5 <^oE|log(l + r,)| 2+5 (6.10) 
where Kq = 2 2+s . Since VjQjVj > 0, Tj (defined in (|6.9jl ) is lower bounded: 

> flV/),^ f 
" " 1 + ^TrZ^-Qj ■ 

Now, since 

< -TrDjQj (-/,) < MttQ,- (-p) < K x = ^=2 sup f- 
n n p n \n 
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and since x i— > is non-decreasing, we have: 



n 3 3 < K 2 = —±- < 1. (6.11) 



1 + ±TrDjQj ~ 1 + Ki 

In particular. Tj > —K-2 > — 1. The function (— 1, oo) 9 a; i— > i2§il±g) i s non-negative, 

A" 2 



non-increasing. Therefore, log ^ +x ) < l°Sil_j£k) f Qr ^ g [— if 2 , oo). Plugging this into {610} 
yields 



E|7j| 2+5 < K K$ +5 'E\T j \ 2+5 

-, 2+5 



= X 3 E|r j | 2+ ' 5 < K 3 E 



y*Q i% - -Tr DjQj 



By lemma 16.21 - (|T]). the right hand side of the last inequality is lower than K± n~( 1+l5 / 2 ) as 
soon as E|Xn| 2+d < 00. This is ensured by A-Q]for S < 6. Therefore, Lyapounov's condition 
(|6.4p holds and Step 2 is proved. 

6.4. Proof of Step 3: Convergence of the normalized sum of conditional variances. 

This section, by far the most involved in this article, is devoted to establish the convergence 
(|6.3p . hence the CLT. In an attempt to guide the reader, we divide it into five stages. Recall 
that z = —p and 

7j = (E j -E j+1 )log(l + r i ) where I, L - Vv , )Q ■ 

In order to apply Theorem l6.1i we shall prove that 0~ 2 Y?j=i ~^ 1 where 2 is given 

by Theorem 13. 11 Since liminf 2 > 0, it is sufficient to establish the following convergence: 



E J+l7? 2 - e 2 . (6.12) 



3=1 

Instead of working with 9„, we shall work with n (introduced in Section[5j see Eq. (|5.ip ) 
and prove: 

n 

X> j+l7 ? - 6 2 0. (6.13) 

3=1 

In the sequel, A will denote a constant whose value may change from line to line but which 
will neither depend on n nor on j < n. 

Here are the main steps of the proof: 
(1) The following convergence holds true: 

n n 

]TE J+l7 2 -£ e j+i (E.T,-) 2 . (6.14) 

3=1 3=1 

This convergence roughly follows from a first order approximation, as we informally 
discuss: Recall that ji = (E,- — E J+ i)log(l + Tj) and that Tj — > by Lemma 
l6.21 (fTj). A first order approximation of log(l + x) yields 7 j « (Ej — E^i)]^. As 
Kj+i(yjQjyj) = —TiDjKj+iQj, one has Kj+iTj = 0, hence 7 j ~ EjTj and one may 
expect Ej + i 7 2 w Ej- + i (E 3 Tj) 2 and even (|6.14[) as we shall demonstrate. 
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(2) Recall that n = E\Xu\ 4 — 2. The following equality holds true 
1 

1 + 



E j+1 (EjTj) 2 = — -1 -2 ( TiDj (E j+ iQj)Dj (E j+1 Qj) 



N 



+ K ^4(E J+1 [Q,] l4 ) 2 \+e 2 j (6.15) 

i=l / 

where 

K 

maxE \e 2 ,j\ < 

for some given _K\ 

A closer look to the right hand side of (|6.15|) yields the following comments: By Lemma [6~3l 
(f2a|) . the denominator (1 + ^TrDjEQ) 2 can be approximated by (1 + ^TiDjT) 2 ; moreover, 
it is possible to prove that [Q^m ~ [T]u (some details are given in the course of the proof of 
step (5) below). Hence, 

N 

^ < 4,(E j+1 [Q J ]uf«-'rmp a . 

n *— ' n J 

i=l 

Therefore it remains to study the asymptotic behaviour of the term -^TiDj(Ej + iQj)Dj(Ej + iQj) 
in order to understand (|6.15p . This is the purpose of step (3) below. 

(3) In order to evaluate ^TrDj(E.j+iQj)Dj(Ej+iQj) for large n, we introduce the ran- 
dom variables 

X < il = ^D^ x Q)D j Q, j<£<n. (6.16) 

Note that, up to rank-one perturbations, EjXj„ is very close to the term of interest. 
We prove here that Xj n satisfies the following equation: 

,., i " -Ti(DiBD m EQ) 1 

xS=- E 7^"V JL ^2x£ n + -TrDeBD J EQ + e 3 , gj , j<£<n, 6.17 

' n m~~j+i (l + ^TrAnEQ) m ' n n 

where B is defined in Section [6. II and where 

maxE|e 3 tj \ < —= . 

t,j<n y/ri 

(4) Recall that we have proved in Section [5] (Lemma 15. ip that the following (determin- 
istic) system: 

n 

m— j + 1 

where a^ m = ^j^yj^^-^yi admits a unique solution. Denote by xf^ n = n (l + ^TrD^T) 2 yf\, 
then (xg*^, j < I < n) readily satisfies the following system: 

p i y ^D t D m T 2 I J<e < n . 
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As one may notice, (|6.17p is a perturbated version of the system above and we shall 
indeed prove that: 

x91 = + £41,- + £42,1 where maxE|e 4 i A < -^= and max |£ 42 ,j | < (5„, (6.18) 
the sequence (<5„) being deterministic with <5„ — > as n — > 00. 
(5) Combining the previous results, we finally prove that 

n 

(EjT,-) 2 - &l 0. (6.19) 

j= i ,wo ° 

This, together with (|6.14j) . yields convergence (|6. 1 3[) and (|6.12p which in turn proves 
(|6.3p , ending the proof of Theorem 13.21 



Proof of (grgp . Recall that ^T^q < ^2 < 1 by (|6~TTj) . In particular, T 3 > -K 2 > -1. 
We first prove that 

E 3 - log(l + T 3 ) = EjTj + enj + £12,3 

where 

In the sequel, we shall often omit subscript j while dealing with the e's. As < K2 < 1, we 
have: 

\k=l / 

fc=2 2 

Therefore, 

Een* < 2lEr^lr,>x : 



(<0 / , ErH, r .|< K2 



(l-^ 2 ) 2 



(b) 2Er1 2Er^ 

< 



( c ) / 2 2 \ / 1 \ 4 (<0 K 

< (^+(r^jE(^--TrD,Q,j < - 

where (a) follows from (a+6) 2 < 2(a 2 +fe 2 ), (b) from the inequality r 2 lr 3 >/f 2 < T 2 Cj^j lr 3 >if 2 , 
(c) from the fact that the denominator of Tj is larger than one, and (d) from Lemma [6T2l pT) 
as Xn has a finite 8 th moment by A-TIJ 

Now, < £12 < EjTjlr^K*- Thus, Eei 2 2 < Er 2 lr 3 >if 2 < K^ 2 ET^lr j>K2 . Lemma 
16. 21 - AT]) yields again: 

„ n .K 
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and (|6.20[) is proved. Similarly, wc can prove 

K 

E j+1 log(l + Tj) = E j+1 Tj + £i3,j with maxE^/ < — . 

Note that since y 3 and F 3 +i are independent, we have E 3+ i(y*Q 3 y 3 ) = ^TrDjMj + iQj which 
yields Ej + iTj = 0. Gathering all the previous estimates, we obtain: 

7j = EjTj + £ll,j + £l2,j - £ 13,j 

= E i r i + £ i4,j , 

where maxj<„ Ee±4j 2 < K n~ 2 by Minkowski's inequality. We therefore have Ej+i(7j) 2 = 
Ey+iQEjTj- +£i4,i) 2 ' Let 

eij = E i+1 ( 7j -) 2 - E 3+1 (E 3 T 3 ) 2 = E i+1 ei 4)i 2 + 2E J - +1 (e 14 , i E^). 

Then 

E| £l)i | < E £l4 / + 2E| ei 4, i E j r i | , 

< E £l4 / + 2(Ee 14 /) 1 / 2 (Er 2 ) 1 / 2 < 



3/2 ' 



where (a) follows from Cauchy-Schwarz inequality and (EjTj) 2 < EjT 2 , and (b) follows 
from LcmmaiUT} which yields ET 2 < Kn' 1 . Finally, we have E"=i E l E j+i (7j) 2 - 
Ej+^EjTj) 2 ] < Err* which implies (pHl)) . □ 



Proof of A6.15p . We have: 



^ = 



1 + iTr^Q 



J Ej (y 3 Q 3 y 3 - ^TiDjQj 



1 + iTr^EQ 



/ J/JQ^i - -Tr£),-Q,- /l 1 
_ E f/j^jtO n j^j -TrD 7 Q 7 - -TrDMQ 

3 \ 1 + ±TtD 3 Q 3 \n 3Wj n 3 H 



Hence, 



E J+1 (E,T,) 2 = - ^i— — E i+1 ( fy*{EjQ 3 )y 3 - ^-TrD^Q, ) + e 21d + e 22 „ 



(1 + iTW^EQ) 2 
(1 + iTr^EQ) 



-L— — E j+1 (y*(E 3 Q 3 )y 3 - -TtD 3 E 3 Q 3 ) + e 2 ,j 
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where 



£21, 



ViQiVi - i TrD,g, f\ 1 
1 + ^TrDjQj \n J ^ J n 3 



1 , 



-2 Ej [ y*Q jVj - -TrDjQj 



E; 



ViQiVi ~ -TrL>,Q, /i i 

n n 



1 + ±TrD 



£ 2,j - 



Ej + l(£ 2 l,j + £22,j) 



(l + lTr^Eg)2 
As ^TiDjQj > 0, standard inequalities yield: 



Ee 2 l,j 



< 



-TrDjQj - —TiDjEQ 



By Lemma [6.21 - (|T|). E(jj*Qjyj — ^TiDjQj) < Kn 2 . Due to the convex inequality (a 
bf < 2 3 (a 4 + 6 4 ), we obtain: 



E -Tr Dj(Qj -EQ) 



E ( — \ v /),:(.) , - EQj) + -Tr Dj(EQj - EQ) 



< KiEl —TtD 3 (Qj - EQj ■ ) + E -TrD 3 (EQj - EQ) 



where the first term of the right hand side is bounded by Kn 2 by (|2cL[) in Lemma 16.31 and 
the second one is bounded by Kn~ i due to the rank-one perturbation inequality [Lemma 
Therefore -.e-j i . / < Kn~ 2 and similar derivations yield E|e22,j| < Kn - '? . Gath- 



ering these two results yields the bound E|s2,j| < Kn 2. Let us now expand the term 
Ej+i (y^EjQjy-j - ^TtDjEjQj) 2 in the right hand side of (jfT2"Tj) . 

Recall that EjQj = E j+1 Qj and that yj = Dj f^Jfv • • > ^ ) T Note also that E J+ i (yjEjQ 
±TrDjE j+ iQj. Then Lemma W% (|2"j) immediatly yields: 



liVi) 



E j+1 ly^QM - —TrDjEjC 

= (Tr^ [E J+l Q 3 )D 3 (Ej+tQ^ + nf^^i <Pi+i[Qj]t 
Equation (|6 . 1 5[) is proved. 



□ 



Proof of (OTP - Recall that x i n = ^TrD t {E j+1 Q)D 3 Q. The outline of the proof of (f67T7]) 
is given by the following set of equations, the x's and e's being introduced as and when 
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required. 

*S =Xi+X 2 -X 3 ( 6 - 22 ) 

Xa = Xa+ea (6.23) 

Xa = X 4 + X 5 +4 (6-24) 

X 5 =X 6 -X 7 + £6-^7 (6.25) 

Gathering the previous equations, we will end up with 

Xt,„ = Xi+X2-X±-X e + X7 + £ where e = -e 3 + £3 - e 6 + e 7 . (6.26) 
Let us first give decomposition (|6.22p and introduce Xii X2 an d X 3 - Recall that B (defined 
in Section [6~Tj) is the N x N diagonal matrix B = diag(frj) where 6j = (p(l + ^TrDiC) 
The following identity yields: 

Q = B + B(B^ 1 - Q~ l )Q 



( -TrAC 



Therefore, 



,U) 



Xi,n = -TrA(E J+1 0)AO 



-TrDgBDjQ + -TrA^diag ( -TrAC* ) (E i+1 Q) DjQ 



1 / ™ 

\m— 1 



j+iymy* m Q DjQ 



A 

= Xi + x 2 - X 3 , 



and (|6.22[) is established. We now turn to decomposition (|6.23p . Identities (|6.ip and 
yield: 

* n — * n — * Q^rnVrnQm _ UmQm _ r * 

VmH — ymHm V m 1 „ — „ " PMmm!/ m Vm ■ 

Using this equation, we have 
X 3 = ^ Tr ^ S ( £ AO 

= ^TrA# ^ E J+1 ([g] mm y m y^Q m )^ AO 
( = } -TrD e B I E c m E j+1 {y m y* m Q m ) J AO 

\m=l / 

-— TrA-B ( E CmEj'+i ([Q]rom \V* m Q m y m - -TrD m EQ J y m y* n Q m ) ) AO 



- X 3 + £3 
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where (a) follows directly from (|6.8p . We are now in position to prove that 



K 

max E | £ 3 1 < —= 

e,j<n Jn 



(6.27) 



Using the fact that \TrAyy* B\ = \y*BAy\ < \\AB\\ \\y\\ 2 together with the norm inequality 
\\AB\\ < \\A\\ \\B\\, we obtain: 



S3 1 < ^^^DjQDtB^Crr&j+lUQ] 



y* m Q m y m — TrL> m E<3 

n 



\ym\nQm\\ , 



(«) a 4 1 ^ 

p 3 n ^ J 



y* m Qmy m — TrAnEQ 

n 



where (a) follows from the fact that \\DjQDgB\\c m < cr^ lax p -3 and [Q] mm ||<9 m || < p~ 2 ■ 
Writing iTr£) m EQ = j^Ti'D m Q + ^Tr D m (KQ — Q) and replacing in the previous inequality, 
we obtain: 



E 



ylnQmVm TrD m EQ 

n 



< 



E 



y* n Q m yrn — TrL> m <3 

n 



2" 


2 






+ 


E 



-TrD m (Q-EQ) 



(3||y„ 



|4\ 2 



where E \y* n Q m ym — ^TrD m <5| < Kn 2 by Lemma ROl -lfTj) combined with Lemma IB~3l 

1 

©, \E\j-TiD m (Q -EQ)| 2 j 2 < ^n-^yLemmaES© andE||y m || 4 < a^ElXu^Nn- 1 ) 
This in particular yields maxf i3 <„ E|£3| < Kn~i and proves (|6.27j) . 
Recall that if m < j, then 



Ej+i(y TO y^Q m ) = E j+1 (y m y^ l )E j+1 (Q m ) = —E j+1 (Q m ). 
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We now turn to Equation (|6.24j) and introduce x 4 , x 5 and £3. 
Xs = ^TcDtB c™%+i (w;Qm)J DjQ , 
= ^TcD t B ( J2 c m D m E J+1 Q m \ D 3 Q 

\m=l / 

+^TrD e B ( a m E J+1 (y m y* m Q m ) ] DjQ , 

= ^TrDzB I J2 £ m D m E J+1 Qj D 3 Q 

\m=l / 

\m=j+l / 

+ -^Tr D e B ^ £ c„,U m E j+1 (Q m - Q)j D 3 Q = Xi + X 5 + 4 . 

and decomposition (|6.24p is introduced. In order to estimate s' 3 , recall that given two square 
matrices R and S, one has iTri^S 1 ! < ||i?||Tr5 for S non-negative and Hermitian. As matrix 
Qm — Q is non- negative and hermitian by (|6.1[) . we obtain: 

141 < \\\DeBDjQ\\ V E(Tr£» m (Q m - Q)) < -™ ax 



(6.28) 



by Lemma 1631 - (f3|. 

We now turn to X5 an d provide decomposition (|6.25[) . Recall that TrAyy*B = y*BAy. 
Combining (|6.ip and (|6.2p . we get Q = Q m - p[Q] mm Q m y m y* m Qm- Plugging this expression 
into the definition \b an d using the fact that y m is measurable with respect to J-j+i (since 
m > j + 1), we obtain: 



X5 = ~^D t B ( ^2 Cm%+1 i.Vmy* m Q m ) ] £>j<2 

.771=7 + ! 



- ^ c m y* m {E j+ iQ m ) DjQ m DiBy m 

m—j + l 

^n[Q]mmym(^j+lQm) DjQmymy^QmDiBy^ 



n 

m=j+l 



In order to understand the forthcoming decomposition, recall that asymptotically y* n A m y m r 
■^TiD m A m as long as y m and A m are independent, and that ^TvD m A m ~ iTr_D m ^4 if A,- 
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is a rank-one perturbation of A. We can now introduce Xe an d X 7 : 



X 5 = - Y —TrD?BD m (E j+1 Q) DjQ 

71 — ' 77, 
L 

n ~2 ^ 

y ^TrAn (Ej+iQ) Dj-Q x -Tr{D m QDtB) + e 6 - e 7 



n * — ' n 

m—j-\-l 

9 n -9 



A 



n ' — ' n 



X 6 - X7 + £ 6 - £7 



where 



m=j+l 



£ V —TrD e BD m (E j+ iQ) Dj 



n * — ' n 

£7 = — y Cm[Q]m m V*n (Ej+lQm) £>jQ 

TO=J + 1 

2 « ~2 1 

V — TrD^ (Ej + iQ)DjQ X —Tr(D m QD^B) 
n t-r 1 n n 

m— j + 1 



It is now a matter of routine to check that: 



K K , 

e 6 < -= and E e 7 < . (6.29) 

\/n Jn 



Let us provide some details. 

Recall that y m is independent from E J+ i(Q m ). To obtain the bound on E|e6 1 , we use 
the facts that E(y* m (E J+1 Q m ) D J Q m D e By m - ±TrD e BD m (E 3+1 Q m ) D Q m ) 2 < Kn" 1 by 
LcmmaOdU), \±TrD e BD m (E j+1 Q m ) Dj (Q m - Q)\ < Kn^ 1 by Lemma , etc. 

In order to prove that E|ey| < Kn~z, we use similar arguments but we also need two 
additional estimates. The control [Q] mm — c m which has already been done while estimating 
£ 3 relies on (|6.8p . The bounded character of E{y^ n A m y m ) 2 where A m is independent of y m 
and of finite spectral norm. This is a by-product of Lemma l6.2I p]l. 
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Wc now put the pieces together and provide Eq. (|6.26|) satisfied by Xij ■ Recall that 

Xi = -TrDiBDjQ 
n 



X2 = 



-TrD e B diag ( -TrAC* ) (E i+1 Q) DjQ 
n \n J 

\m<j J 

X 6 = ^TiD^bI J2 c m D m \{E J+1 Q)D 3 Q 

\m=j+l J 
2 n ~2 i 

P_ V- — TrL> m (E j+1 Q) DjQ x -Tr(D m QD e B) 

m—j-\-l 

n ^ (1 + ±TrD m EQ) 2 Xm ' j ' 

As n ELi d rnD m = diag (±TtD x C , ±TrD N C^ , we have x 2 - Xi - X 6 = 0, and 
becomes 

m 1 1 A -T^(DfBD m Q) 

m=j + l y n m ^6) 

where E|e| < Kn~^ thanks to Inequalities ()6.27|) . ()6.28|) . and (|6.29|) . Small adjustments 
need to be done in order to obtain (|6.17p . Now replace —TtD(_BD p Q by -^TrDgBD p EQ (use 
Lemma IOI (|2^|) ) . The new error term £z,tj still satisfies maxjj<„ E|e 3i ^ )3 -| < Kn~*. Eq. 
(|5T7|) is proved. □ 

Proof of \6.18\ ). Recall that xt,j an d £3,£j have been introduced above. 

Following the matrix framework introduced to express the system satisfied by the y's (ma- 
trices An, and A n '^), we introduce matrix G n = Aj t = {ge m )i m =u its (n — j + x (n — 
j + principal submatrix Gn = {ge, m )e m =j an< ^ the matrix Gn which differs from matrix 

Gn by its first column, equal to zero. Denoting by 5^ = ^TtDiDjT 2 ; j < t < n), we 
have: 

x W=G°MxW + 6W . 

Define here the (n-j + l)xl vector £3 = (£3,^; j < I < n) and the (n — j + 1) x 1 
vectors: 

x (j) = (xS; i < t < «) . 

Define now the (n — j + 1) X (n — j + 1) matrix: 
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and G°'U) which is egal to G^i exept for its first column equal to zero. With these notations, 
Eq. (|6.17p . valid for for j < I < n, can take the following matrix form: 

We will heavily rely on the following: 



lim sup 



j _ qO,U) 



< OO . 



which can be proved as in Lemma T5. 21 - (|4|) and Lemma T5. 51 We drop superscript 0, V) in the 
equation below for the sake of readability. 

X = Gx + $ + £3 

«*• X = G X + S + e 3 + (G-G)x + (S-S) , 

& X = (I- G) -1 * + (/ - G)~ 1 e 3 + (J - G)-\G - G) X + (I - G)' 1 ^ - S) , 

O x = x + (I-G)- 1 e 3 + (I-G)- 1 (G-G)x + {I-G)- 1 (S-S) (6.30) 

The first component of the previous equation writes: 



4 J) + [(/-G 'W)- 1 e 3 ]i 

+ [(/ - G r °'0'))- 1 (G '^ - G a - ij) )x + G ' U) y\S - 6)] j 



A (j) 

= a;. + £41 j + £ 42 j 



Due to Lemma l5?2l -(j4| which applies to G '^ and to the fact that max<j<„ Eles,^- 1 < Kn 2 . 
we have: 



n-j+l 



The second error term £42^' is the sum of the following terms: 

£42 tj = [(I-G '^)- 1 ^'^ - G°' U) )x]i + [(I - G ^)- 1 ^ - 6))t 

Let us first prove that [(I — G '^) -1 (G '^ — G°'^)x}i is dominated by a sequence indepen- 
dent of j that converges to zero as n — > 00. The mere definition of Xitn ( see Q6.16[) ) yields 
llx^'lloo < (-^ cr 4ax)( n P 2 )~ 1 j where || • ||oo stands for the ^co-norm. Hence 

- G°.U))-i(G°.«) - G°' U) )x]i\ 

< (7-G '^)- 1 (G°<«-G°<W) T \\x\\oo<K (G°<«-G°<«) T 

OO OO 

Let us prove that 



-> 



00 n— >oo 



(6.31) 
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uniformly in j. To this end, let us evaluate the {£, m) -element of matrix G '^' — G '^' 
(to > j): 

n\[G°^-G°^] em \ - 



±TrD e D m T 2 



< 



< 



(1 + ±Tr£ m EQ) 2 (1 + iTrAnTT 
(1 + -TiD m T) 2 -TiD e BD m EQ - (1 + -TrD m EQ) 2 -TrD e D m T 2 



(1 + — TrD m T) 2 —TrD e BD m (EQ - T) 
n n 

(1 + — Tr D m T) 2 — Tr D{TD m (B - T) 
n n 



(1 



-TrD m T) 2 - (1 + -TvD m EQ) 2 
n n 



1 



TvDfDrnT 2 



(6.32) 



The first term of the right hand side of (|6.32[) satisfies: 



(1 + —TrD m T) 2 —TrD e BD m (EQ - T) 
n n 



< 1 



-Tr£/(EQ - T) , 



where U is the N x N diagonal matrix U = cr^ax P 1( hag (sign {E[Q]u — ti) ,1 < i < N). By 



Lemma 16.31 (|2a[) . the right hand side of this inequality converges to zero as n — ► oo. 

The second and third terms of the right hand side of (|6.32|) can be handled similarly with 
the help of Lcmma l6.3l and one can prove that elements of n(G '^ — G°'W) are dominated by 
a sequence independent of j that converges to zero. This implies that (G 0, ^- ) — G°'^) T 

oo 

converges to zero uniformly in j and (|6.31[) is proved. As a consequence, [(J— G '*-^) -1 (G 0,( ^ — 
G '")^]! is dominated by a sequence independent of j that converges to zero. The other 
term [(/ — G°'W) _1 (5 — S)]\ in the expression of £42,j is handled similarly. Eq. (|6.18|) is 
proved. □ 



Proof of UTW\) . We rewrite Equation (|6~15)) as E j+1 {EjTj 
1 1 



n 2 h 



^TrDjEQ) 
1 



: TJi j + KT]2,j + £2,j with 

^TrDjiEj^Q^D^E^Qj) , 



N 



^l + iTr^EQ) 2 ^ ^ 



V n -—^ and 2j=i — W n where V„ and W„ are 



and we prove that Y]j—i 

defined in Section [5] To prove the first assertion, we first notice that 
TtDjfEj+tQADj&j-nQj) = E J+1 (1iD j (E j+1 Q J )D J Qj) = E J - + i(TM> i (E i+ iQ)D i Q) + e 



with lei < 2at 



by Lemma 



Therefore, 



'/i, 



iTr_DjE<5)' 



It remains to control the difference (l + iTrl^EQ) 2 - (l + ±TrDjT) 2 , to plug (f67T8]) 
and one easily obtains Xw=i ?7i,j — V„ -^-> 0. 



3S 
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We now sketch the proof of X^j=i V%,j — ^« — * 0- m f)6 . 2 [) . satisfies 



— (z(l+£jQl£i )) where ^ is the row without element j, and Q\ = (Y- Y? + pl„_i) 
where YP is matrix Y without row i and column j. Using this identity and Lemmas 16. 21 (TT1) 
and 16. 3[ we can show that [Qj]« is approximated by £j, which is key to prove Y^j=i V2,j ~ 

W n ^0. □ 

7. Proof of Theorem 13.31 

We first provide an expression of the bias that involves the Stieltjes transforms -^TrQ 
and jjTrT. By writing log det(Y n Y* + pi N ) = Nlogp + \ogdet(^Y n Y* + I N ) and by taking 
the derivative of logdet(-Y^i^ + Zjy) with respect to p, we obtain 

logdet(F„y„* + P I N ) = TV log p + N (J- - 1ttQ(-o;)J du . 

Since ^TrQ(z) € <S(R+), we have £ - ^-TrQ(-w) > for lu > 0. In fact, recall that 
IIQ( — w )| < by Proposition ^. 21 Therefore, by Fubini's Theorem, 

Elogdet(F„r„* + P I N ) = Nlogp + N J (J- - ^TrEQ(-w)^ dw . 

Similarly, 

NV n (p) = N J log(A + p)vr„ (dX) =N log P + N J - -^TrT(-w)^ dw. 
Hence the bias term is given by: 

B n (p) = Elogdet(r„y„* + pI N ) - W n (p) = / Tr (T(-w) - EQ(-w)) cJw . (7.1) 



In Appendix IB1 we prove that: 

Tr(T -Q) = Tr(f - Q) . (7.2) 
Therefore, we can also write the bias as: 

B n (p) = [ Tr (f(-w) - EQ(-w)) cJw . (7.3) 
Jp 

For technical reasons (and in order to rely on results of Section [5|) , we use representation 
(|7.3[) of the bias instead of (|7.1[) . The proof of Theorem 13.31 will rely on the study of the 
asymptotic behaviour of the integrand in the right hand side of this equation. 

As a by-product of Section [5j the existence and uniqueness of the solution of the system 
of equations (|3.1[) is straightforward. Indeed, define the n x 1 vectors w and p as 

W = (Wj. n ; 1 < j < n) , 
P = [Pj,n\ l <3<n) . 
Then the system (|3.1I) can be written in a matrix form as 

w = Aw +p . (7.4) 
Since (I — A) is invertible for n large enough, this proves Theorem 13. 31 - (fTj). 
The rest of the proof will be carried out into four steps: 
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(1) We first introduce a perturbated version of the system (|7.4p . For the reader's con- 
venience, we recall the following notations: 
1 



U 



U 



h = 



^(1 + iTi'AT 
1 

1 + iTr AEQ 
1 



^(1 + iTr^T) 
1 



(l + ilM^EQ)' 

1 

[1 + lTrDjCy 



where z is equal to — u> with uj > 0. Write the integrand in (|7.3[) as 
1 n 

Tr (f{-u) -Wii-wj) = -^cp^u) with tp^u) = »&(-«<;) - U[Q(-u)]jj) ■ (7.5) 

71 3=1 

Let tp^ (uj) = n(bj(—uj) — K[Q(—Lj)]jj) and define the n x 1 vectors (p and ip and 
the n x n matrix A as 

¥> = (Vj-; 1 <j < n) , 

if) = (ip u) ;l < j < n) , 



A : 

We first prove that 

(2) We prove that 
^,00 



±TrDjD m CT 



(l + lTr^T)(l + lTr^C) 
(p = A(p + ip . 



j,m=l 



(N / n \ _ iV > 

£ E 4^ E a im E lQi\lim - ^ E 4 K ^1: 
i—l \ rn—1 / i—1 / 



t (3) 



(7.6) 



(7.7) 



with |e^| < Kn where if is a constant that does not depend on n nor on j 
(but may depend on u>). 

(3) Matrix A readily approximates A and vector if) approximates p for large n by Step 
[2] Therefore, by inspecting Equations (|7.4[) and (|7.6[) . one may expect <p to be close 
to w. We prove here that 



ll^-HU »0. (7.8) 

n—*oQ,N/n — 'c 

(4) Let p n {w) = i jy; =1 w j>n (u). Eq. yields A £» =1 - /?„(cj) . In order 

to prove (|3.5p . it remains to integrate and to provide a Dominated Convergence 
Theorem argument. To this end, we shall prove that: 

l&MI < ^ (7-9) 
for n large enough. This will establish (|3.4p . We will also prove that 



E^(° 



< 



A"' 



(7.10) 
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for oj G [p, +00), where K' does not to depend on n nor on to. This will yield (|3.5[) 
and the proof of Theorem 13.31 will be completed. 



7.1. Proof of step QJ Equation (^6]> . Recall that W = n(6,- - E[Q]jj)- Usin S thcsc 
expressions, we have for 1 < j < n: 



(Pj = n{t ] -b J )+tl} (j) = nbjtj (bj 1 - ij 1 ) + ipW 
= LubjtjTrDj (C — T) + (j) 

N 

GO 



LO 

i=l 



to 



N n 



n 

i—l m— 1 



n 



n 

m—l 



which yields Eq. 



7.2. Proof of step [2) Expression of i/j^*'. We shall develop ijyi' as 



$W = X + 2 - 3 (7.H) 

^i=^4 + ei (7.12) 

^ 2 = --0 5 +-0 6 (7.13) 

5 =0 7 + £ 5 (7.14) 

6 = 8 + £ 6 (7.15) 

03 =0 9 + £3 (7.16) 



where the fc 's and the e^'s will be introduced when required. We shall furthermore prove 
that |£ fe | < Kn- 1 ' 2 for k = 1,3,5,6. This will yield 



W = 4 - 7 + 8 - 9 + e« with | £ 0')| = | ei - £3 - e5+e6 | < . (7.17) 
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Let us begin with decomposition (|7.11[) : 

tfCi) = nbjE ([Q]jj ([Q]jj~ — ^J 1 )) 

( = } nubjE ( [Qjjj (y* Q 3 y 3 - - TrDj C 



= nuibjCjE (vjQjVj Tr DjC 



bfc 3 E [Q]jj y*Q 3 y 3 - -TrDjEQ y*Q jVj - -TrDjC 



= LobjCjTrDjE {Q 3 - Q) + ubjCjTrDj (EQ - C) 



-nu'bjCjE \\Q\jj [y*Q jyj - -TrDjEQ j \y*Q 3 y 3 - ~TrD 3 C 

= V>1 + ^2 ~ ^3 

where (a) follows from (|6.2[) and the definition of b 3 , (b) follows from identity (|6.8[) , and (c) 
follows from the following equality: 

E (vjQjVj - ~ Tr£) J C ') = ^ Tr£) J ( E( 5j ~ C ) = ^ Tr ^' E ^ - ^) + ^ TrZ ^ ( E<9 ~ G "> ' 

Eq. (|7HT) is established. 

We now turn to the decomposition (|7.12p . Combining (|6.1[) and (|6.2j) . we obtain Q = 
Qj - Lu[Q] jj Q j y j y*Q 3l hence Vi = ^ 2 b 3 c 3 E {^Q] 33 y*Q 3 D 3 Q 3 y 3 ^j . Using identity (|6^8|) and 
the fact that E(y*Q 3 D 3 Q 3 y 3 ) = ±E(Tr DjQjDjQj), we obtain: 



*3 

r 
n 



Vi = ^6 3 - £j E (TrD 3 D s Q 3 ) 



•%<fE f [Q]# [y*Q 3 y 3 - ^TrDjEQ^ {y*Q 3 D 3 Q 3 y 3 )^j 



A 



We have: 



|ei| < —E (y*Q 3 D 3 Q 3 y 3 \e n + £12 + £i 3 |) , 



with e u = y*Q 3 y 3 - j- l TrD 3 Q j , e 12 = ±TrD 3 (Q 3 - EQ,), and e 13 = ^TvD 3 E (Q 3 - Q). 
By Lemmas 0)3,0-0, and0-©, wc have E|e n | 2 < KrT 1 , E\e 12 \ 2 < Kn~ 2 , and 
I £■ 13 1 2 < Kn~ 2 respectively. By Cauchy-Schwarz inequality, we therefore have: 



K(E{ V *Q 3 D 3 Q 3 y 3 fy- K< 
|£l1 " ^ ~JK 



and (I7.12|l is established. 
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We now establish decomposition (|7. 13|) : 
■02 = ubjCjTrDj (EQ - C) 

N 
i=l 

N 



-uj 2 



(a) 



-C 2 



i=i ^ ^ 

Sigj E 4 c i E (fitter) - -i^aeq ) 



i=l 



-^5 + ^6 

where (a) follows from (|6.T[) . Equation (|7.13|) is established. 

Let us now turn to decomposition (|7.14[) . We have 5 = - b ^ Cj J2i=i ofj^TtDiE, f*9i — Q 
By similar arguments as those used for we have: 



^ = -^E4^ E ( Tr ^^ 

w 6,'C 



1=1 

JV 



i=l 



A 

= 07 + £ 5 



where |es| < Kn s and (|7.14p is established. 
Turning to (|7.15|) . we have 

w / / 1 \ 2N 



06 = w 3 5 3 -c,- £ 4c 2 E i[Q]u (&Qi& - yTvD^Q 
N / 1 



^JCj E4 C ? E - ^TrAEQ 



A 



06 + £ei, (7.18) 
using again (|6.7[) . The term £6i satisfies: 

1 ^ 3 
|S61 1 < — y^g t 2 7 E | £611,i + £612, i + £613, d 



2=1 



q N 

< ~Y. a l ( E l e 6H ) i| 3 + E |£612 > i| 3 + |e613 ) i| 



LO 

i=l 
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where e 6 i M = -^AQi, e 612 , ? = ^TrA [Qi - EQi ) , and e 613> i = ^AE [Qi - Q 



ByLemmaiS(II]),E|e 6 ii, ? ;| 3 < Kn^' 2 . By Lemma[0]-([2d]), E |e 612 ,i| 3 < (E|e 6 i 2 , 
Kn~ 3 / 2 . By Lemma , |£ 6 rm| 3 < ifn~ 3 , hence 

K 

£6,1 < -7= • 



3/4 



< 



We now handle the term ip' e in (|T. 18[) . We have: 

N ( 1 

i/>' 6 = ujHfcj cr 2 c 3 E I - -TiDiQi + e 6 i2,< + £ 6 i3, 

i=l ^ n 

iV 



2 c 3 J 



■08 + £ 62 , 



where 



£62 = W 



2 c 3 



T^DiQi +e 62 

n 



1 



(e 6 i2,i + e 8 i3,i) + 2E fiQiC* - -TrAQi (e 8 i2,< + eeia, 



Using Lemmas 16.21 - (fT|) and 16.31 it is easy to show that 

K 

£62 < -f= • 



Furthermore, the terms E ( ) 2 in the expression of i/> 8 has a more explicit form. Indeed, 
applying Lemma |6~2"1 P|) yields: 

= E a M E ( Tr ^QiDiQi) + k E ([Qi]^) • 

i=l \ m=l / 

Decomposition (|7.15p is established with eg = £gi + £62- 

It remains to give decomposition (|7.16p . Using (|6.8[) , we have 

^ 3 = no. 2 6 J 5 J E ([Q)jj (y*Qm - -TtDjEq) (y^QjVj - -Tr DjC 



nw%c 2 E ( ( HjQ,,/, - -TrDjEQ I ( y*Q m - -TrD i C' 



-nw 3 6j-c 2 E 



J j j 



y*Q jVj - -TrDjEQ tfJQjty - -TiDjC 



A 



^3 + £ 31 

The term £31 satisfies 



|£ 31 | < — E 



VjQjVj TtD jQj + £ 311 + £ 312 



VjQjVj — TTD jQj + £ 3n + £312 + £3 



13 
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with e 3U = ^Dj(Qj - EQi), £312 - ^TrDjEiQj - Q), and e 3 i 3 = ^TrDj(EQ - C). 
The terms £311, £312, and y*Qjijj — ^TiDjQj can be handled by Lemmas l6.2K [Tj) and 16.3 
The term £313 coincides with ■0 2 ( Tlw ^iCj) _1 - The derivations made on 2 above (decompo- 
sitions Q7.13H7.15]) ) show that \ip 2 ( w bjCj)~ 1 \ < K therefore | £r 313 1 < KrT 1 . 

Using these results, we obtain after some standard manipulations: 

K 

631 ^ -77, ■ 

The term 3 can be written as: 

3 = nw 2 6j-cfE ( (y*Q jyj - ^TrDjQj + e 3n + e 3U 

y*iQjVj - -TzDjQj + £311 + £312 + £313 

n 

nuj%£ 2 M ( y*Q. m - -TrDjQj) + £ 32 



3 \ »i 

A 



77 



■09 + £32 

with |e 3 2| < KnT 1 / 2 . Similarly to 8 , we can develop 9 to obtain 



^~ h A ( ^ 4 , - 

09 = ^DjQjDjQj) + «E 4 E {[Q& 



(7.19) 



i=l 



Decomposition f|T. 16[) is established with £3 = £31 + £32. 

We now put the pieces together and provide Eq. (|7.17p satisfied by 0^ . We recall that 



t/> 4 = — E (TrDjQjDjQj) , 

uj 3 b ~ N 

^7 = i ^p I ^4c 3 E(TrAQ l AQ l ) , 
i=\ 

3l ~ N ( 

09 = — 1(1W¥?^Q;) + "£4 E ([Q 



2 

mm 



jV 



When computing the right hand side of (|7.17[) . all terms of the form KTiDjQjDjQj and 
ETrDiQiDiQi cancel out and we end up with Equation (|7.7p . Step is established. 
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7.3. Proof of step [3j \\ip-w\ 



0. In order to prove 



(A-Af 
limsup (I 



oo n— »OQ 

1 III 



AY 



< oo, 



, we need the following facts: 

(7.20) 
(7.21) 



I — A is invertible for n large enough, 



lim sup 



(i -Ay 



< oo 



(7.22) 
(7.23) 



The proof of (|7.20p is close to the proof of (|6.31[) above and is therefore omitted. The bound 
([7721]) follows from Lemma [572f-([3)l. We now prove ([7722]) and ([7723)1 . Recall that by Lemma 
15.51 there exist two vectors u n = (ui, n ) >- and v n = {v^ n ) >~ such that u n = Au n + v n , 
sup n ||u n ||oo < oo and liminf n min£(f£_ n ) > 0. Matrix A satisfies the equation u n — Au n -\- v n 



with Vr, 



(A — A)u n . Combining (|7.20p with inequalities sup,. 



< oo 



and lim inf „ (min£ V£ n ) > 0, we have liminf„(min^ vg n ) > 0. Therefore, Lemma 15.21 applies 
to matrix A for n large enough which implies l|7.22[) and (|7.23D . 



We are now in position to prove \\ip — w 
obtain: 



-> 0. Working out Eq. dLSj and ([774]) . we 

<p = w + (I - A)" x (i - A)<p + A)-\iP - p), 



hence 



lv-Hloo< - AyHl (A-A) || ¥ ,|| 00 + |||(7-^)-i||| iiv-pII 



Thanks to (|7.22j) . we have ip = (I — A) 1 xp for n large enough. One can check from (|7.7[) 
that sup n Halloo < oo. Therefore, by ()T. 23[) . we have sup n Halloo 

< oo. Using ([7720]) and 

([7721]), we then have III (I- A)" 1 III {A-A) WtpW^ -> 0. 

00 oo 

It remains to prove that — Plloc —> 0. In Step 3, it has been established that if) is a 
perturbated version of p as defined in (|3.2p in the sense of Eq. (|7.7[) . Using the arguments 
developed in the course of the proof of (|6.18[) . it is a matter of routine to check p||oo — * 0. 
Details arc omitted. Hence 



Consequently, \\<p — w\ 



WKi-Armjw-pw 

and Step [3] is proved. 



0. 



7.4. Proof of step [4j Dominated Convergence. In this section, constant K' does not 
depend on n neither on uj but its value is allowed to change from line to line. We first prove 
(fT9]) . Wc have 

i/su < Mi-. ^ lll^-^r'IILiipiu 

by ([773]) . By inspecting ([372]) one obtains UpU^ < |/c|(^/n)(^p + ^p) < K'uj" 3 . Wc need 
now to bound |||(7 — ) 1 1 1 1 in terms ofw€ [p, oo). Lemma [5721 - ([73]) yield: 

|||(7- A)- 1 III < max ^^.") 
00 ~ min^(^,„) 
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where u n = (u£ n ) and v n = (vi n ) are the vectors given in the statement of Lemma 15.51 We 
now inspect the expressions of ue n and vi n . Eq. (|5.4[) yields: 

min(u £! „) > — ^ -2 min-^TrDj, 

and max^(it^ ra ) < (Na^ aax )(nui 2 )~ 1 by (|5.6|) . Gathering all these estimates, we obtain |/? n | < 
if co~ 3 . and Inequality (|7.9[) is proved. 



We now prove (|7.10[) . We have 
1 " 

i-5>ji<iMioo< c-^r 1 ii^ii< 



(7.24) 



by (|7.6p and (|7.22p . We know that the right hand side is bounded asn^oo. However, not 
much is known about the behaviour of the bound with respect to lo. Using Inequality (|7.24|) 



and relying on the derivations that lead to (|7.6H7.7p . one can prove that (I — A) 1 

CXJ 

H^Hoo, and therefore ||v|l°o are bounded on the compact subsets of [p,+oo). Therefore, 
in order to establish (|7. 10|) . it is sufficient to prove that ||y||oo is bounded by K' ui~ 2 near 
infinity. To this end, we develop |<p-(o;)| as follows: 



\vM\ 



lit ; 



e {[Q\n {[Q]Ji tj 1 



nuitj 



< wijT&lQ}^ \TiDjE (Q - T)\ + nujij 

(o) _ „ 

< utjElQljj \TrDjE (Q - T)\ + ncotjCj 



E ( [Q]jj ( VjQiVi - -TrDjEQ 



E y*Qm - -TrDjEQ 



+nuj 2 tjCj 



E [Q]jj V*Qm - -TrDjEQ 



< utjE[Q]jj \TtD 3 E (Q - T)\ + LoijCj \E(TvD J (Q j - Q))| 



+2nw 2 i j c j E ( [Q] 3J [y*Qjyj - ^TrDjEQj 

. UJ 2 tiC 



+2—±l.E[Q] Jj {TrDjE(Qj - Q)) , 

where (a) follows from (|6.8p and (b), from the fact that 

{y*Q jVj - ^-TrDjEQ) 2 < 2(y*Q jyj - ^-TrDjEQj) 2 + 2(^-TrD j E(Q j - Q)) 2 . 

Let ot(u>) = fimaxi<,<]v \U — E[Q]u\. Using Lcmma l6.3l (|3j). we obtain from the last inequal- 
ity 



vMlloc < ^s*a(w) + G ^f- + — E (y*QjVj - -TrDjEQ 



2ai 



As in (|7.19p . we have 



E 



JV 



ViQiVi ~ — TrDjEQ j = - E(TrDjQjDjQj) + k £ ^[QiM < 



Na^{l + \n\) 
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Therefore, 

\\(p(u)\\ 0O <^SLa(u) + ^ (7.25) 

2 , 

for lo £ [p,+oo). A similar derivation yields a(u) < ||y(kj)||oc + Plugging this 
inequality into (|7.25p . we obtain 

(i-^/^ii^miu < ^, 

hence ^(a^Hoo < K'lo~ 2 for co large enough. 

We have proved that ^(w)!^ is bounded on compact subsets of [p, oo), and furthermore, 
that (|7.10p is true for ui large enough. Therefore, (|7.10p holds for every ui E [p, +oo). StcpQ] 
is proved, and so is Theorem 13.31 



Appendix A. Proof of Lemma 16.31 
Proof of Lemma l6.3K}Tj) . Straightforward. 
Proof of Lemma ROl - © . 

Proof of H2a\). From (TTJ, Lemmas 6.1 and 6.6], we get 

-TiU{Q(-p)-T(-p)) >0 a.s. 

Tl n— *oo 

Now since 

2||£/|| 



-TvU(Q(-p)-T(-p)) 
n 



<ll^ll(IIQ(-p)ll + l|r(-p)ll)< 



p 

the Dominated Convergence Theorem yields the first part of (|2a[) . The second part is proved 
similarly. □ 

Proof of IW\). Recall from Theorem 12. 3K 1) and from the mere definitions of T and B that 
matrices T(z) and B(z) can be written as 

f-V Dj ] and B= -zI|-V =-i D, 

n Z^l + L TrDjT J J y n^l + iTr^EQ J 

We therefore have 

-TrU(B(-p) - T(-p)) = -TrUBTiT- 1 - B" 1 ) 
n n 



1 / <A -TrDjfEQ - T) 
— Tr UBT^p- ^ V, , -£> 

n2 ^ ^(l + iTr^T)(l + lTr^Eg) 

^ N n 



i=l j=l 



[WubiUo-fi 1 

with a;™ = , L - L — -\ r—TrDJEQ - T). It can be easily checked that 

13 (l + iTr^-T)(l + iTr^-EQ)n jV V 7 ' 

< 2sup„(||C/||)cr4ax//° 3 - Furthermore, x" — >„ for every i,j by (|2"aj) . It remains 
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to apply the Dominated Convergence Theorem to the integral with respect to Lcbcsgue 
measure on [0, l] 2 of the staircase function f n (x,y) defined as f n (i/N,j/n) = xfj to deduce 
that ^TtU(B - T) -> 0. This ends the proof of ([2b]). □ 



In the sequel, K is a constant whose value might change from line to line but which 
remains independent of n. 



Proof of [2qj . We have 

(a) 



W(Q-EQ) l = J ^2(Ej-E j+1 )TrUQ 



3=1 



(b) 



^(E j -E j+1 )TrC7(Q-Q i ) 

3=1 



H -£<%-%«>f^gf ft X>, (a..) 

3=1 L + yjViVj j=1 

where (a) follows from the fact that EiTr[/Q = TrZ/Q and E n+1 Tr UQ = A vl (). (b) 
follows from the fact that EjTiUQj = Ej + iTr UQj since Qj does not depend on yj and (c) 
follows from (|6.1[) and from the fact that Tr QjVjVjQjU = VjQjUQjyj. 
Now, one can easily check that J2j=i x j ( = Tr?7(Q — EQ)) is the sum of a martingale differ- 
ence sequence with respect to the increasing filtration !F n , . . ., T\ since E^Xj = for k > j. 
Therefore, 



E{TtU(Q-EQ)) 2 = J2 

3=1 

Write where: 

VjQjUQjyj 



Xj,i = -(Ej-Ej+i) 
£j,2 = -(Ej-Ej+i) 



1 ' - 

VjQjUQjyj VjQjUQjyj 



1 + tiQiVi 1 + ^DjQj 



Using the fact that yj and Fj+i are independent, and the fact that Qj does not depend on 
yj, one easily obtains: 



VjQjUQjyj \ l / Q,[/Q 



J+1 U + ±TrD j Q j ) n J 3+1 Vl + iTr^Q, 
Thus Xj^i and write: 

* \1+ ^TiDjQjJ n 3 3+1 {l + ^TrDjQ, 

m w \ VjQjUQjVj ( * n i Tnn 

Xj2 = (E.i — E.; + i) ; ; — VjQiVj llUMj 

= (Ej - E j+1 )xj t3 . 
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Since matrix ||£)jEj ( y 



Q 3 UQi 



L+iTrDjQj 



< K, Lemma 16.21 - (fTj) in conjunction with Assumption 



A-Q] yield Ex\ < Kn 1 . Furthermore, we have: 



X 3,3\ < 



VjQjUQjVj VjQjVj - -T^DjQj 



since y*Qji/j > and -^TrDjQj > 0. Cauchy-Schwarz inequality yields: 



E^ 3 < {E{v]QjUQ ]yj T) ME y*Q jVj - -Tr DjQj 



which in turn yields Ex 2 3 < ^ since 

E^QjUQm) 4 < K and E [y]Q jVj - ^D.Q^j < -J , (A.2) 

where the first inequality in (|A.2|) follows from < VjQjUQjyj < \\QjUQj\\ \\yj\\ 2 and from 
Assumption A-Q] and the second from Assumption A-Q] and Lemma l(T2T (fTj) . 

We are now in position to conclude. 

Ex\ 2 = E((E j -E J - + i)x i>3 ) 2 < 2E ({E jXj , 3 ) 2 + (E J+lXj , 3 ) 2 ) 

< 2E(E,-4 3 +E i+1 4 3 ) = 4Es£ 3 , 
where (a) follows from Jensen's inequality. Now, 

Ex] = E(xj,i + x ji2 ) 2 < ((Ea£i)* + (Ea$ 2 )*) < ^ 
and E(Tr£/(Q - EQ)) 2 = £" =1 Ex 2 < K. Inequality © is proved. □ 



Proof of k2d\) . We rely again on the decomposition (jA.lj) and follow the lines of the compu- 
tations in ([3], page 580): 

TrU(Q-EQ) = -±(E j -E j+1 ) y J^ jV 

j—l 1 + Vi 



Thus, 



iVi 



(«) if 



ViQjVj 



K ( , ,y*jQjUQjy., 
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where (a) follows from Burkholder's inequality and (b) from the convexity inequality (J^iLi cii) 
n y^Lj Recall now that y*QjUj > and ||Qj(— < Standard computations yield: 

y 1 t Uj Wj yj ) P 

which is uniformly bounded by Assumptions A-Q] and A-[2] Therefore, (|2d[) is proved. □ 



Proof of Lemma !6.3Kj3|) . Developing the difference Q — Qj with the help of 
obtain: 



|TrM(Q - Qj 



TtM 



Qj H :!/'/> j 



\yjQjMQjyA 



< \\M\ 



\QjVj\ 



Consider a spectral representation of Y^Y^* , i.e., Y' J Y'i* = Y^,i=i XiCie*. We have 

,*„, . |2 



AT 



N I * |2 

^mt^TTT^ and y*Q jVj = Y. 



hence the result. Inequality ^ is proved. 



we 



Appendix B. Proof of Formula (|7.2|) 

Recalling that Q(z) = (YY* — zl^)^ 1 and Q(z) = (Y*Y — 2J„) _1 , it is easy to show that 
Tr(Q) - Tr(Q) = (n- A)/z. We shall show now that Tr(T) - Tr(f) = (n - iV)/z. Formula 
(|7.2|) is obtained by combining these two equations. 
Equations (|2.2j) in the statement of Lemma 12.41 can be rewritten as 

t ™ 1 t N 1 

** + ^ E 4*i = - for 1 < * < ^ *i + i E 4*< = - for 1 < j < n - 

3=1 i=l 
By summing the first iV equations over i and the next n equations over j and by eliminating 
the term i 2j=i a ijtitj> we obtain yj\ £j — y\ tj = [n — N)/z, which is the desired 

result. Equation (|7.2p is proved. 
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